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Abstract 

Some of the theoretical challenges posed by the general relativistic 
description of binary systems of compact objects (neutron stars or black 
holes) are reviewed. We recall the various ways one can use the theory of 
the motion, and of the timing, of binary pulsars to test the strong-field 
and/or radiative aspects of General Relativity. Recent advances in the 
theory of the motion and radiation of binary black holes are discussed. 
One emphasizes the usefulness of the Effective One Body approach in 
providing a quasi-analytical description of the waveform emitted by coa- 
lescing binary black holes. 

1 Introduction 

The discovery, in our Galaxy, of binary systems comprising gravitationally con- 
densed objects (neutron stars or black holes) has opened up both new exper- 
imental opportunities, and new theoretical challenges. Here we shall focus on 
the theoretical challenges posed, for certain binary systems, by the necessity of 
getting a very accurate, general relativistic description of binary systems made 
of condensed objects. We have particularly in mind two different physical situ- 
ations. 

On the one hand, the discovery of binary pulsars in 1974 pQ has given us the 
challenge of developing a theory of the relativistic motion of two compact objects 
which is accurate enough to match the remarkable precision of the observational 
data. Indeed, the very high stability of 'pulsar clocks' has made it possible to 
monitor the orbital dynamics of binary pulsars down to a precision allowing one 
to measure secular effects linked to very small (~ (v/c) 4 and ~ (v/c) 5 ) terms in 
the orbital equations of motion, as well as periodic effects linked to 0({v/c) 2 ) 
terms. An additional challenge is that these small 'post-Newtonian-type' terms 
in the equations of motion must be cleanly disentangled from the numerically 
much larger self-gravity effects. Indeed, though both types of effects can be 
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formally expanded in a post-Newtonian (PN) expansion in powers of 1/c 2 , the 
self-gravity effects contain powers of 7,; = GM/(c 2 R), where R denotes the 
radius of one of the compact objects, while the orbital effects should gather 
only the terms containing powers of j e = GM/{c 2 D) ~ (tWbitai/c) 2 , where D 
denotes the typical distance separating the two objects. [The subscript i in 7$ 
refers to 'internal', by contrast to the subscript e which refers to 'external'.] 
One can always consider 7 e as a small parameter (for instance j e < 10~ 5 in 
the currently observed binary pulsars), while, because R -C D, one will have 
7i 3> 7e- For a usual star, or even a white dwarf, 7, is quite small and can 
be used as an expansion parameter. By contrast, for a compact object (here 
defined as a neutron star or a black hole) 7$ will be of order unity, and cannot, a 
priori, be meaningfully used as an expansion parameter. We shall review below 
some of the methods which have been used not only to show how one can, in 
principle, disentangle the 'orbital' (j e ) expansion, from the 'self-gravity' (7*) 
one, but also to compute the 7 e expansion to the very high accuracy needed to 
discuss observational data. 

The second physical situation which yields an even bigger theoretical chal- 
lenge is the forthcoming detection of gravitational wave signals, in large inter- 
ferometers (LIGO, VIRGO, GEO600, LISA,...). Indeed, one of the premier 
sources that one hopes to detect in LIGO/ VIRGO is a coalescing binary black 
hole. The observationally important signal for these sources is generated during 
the last few orbits leading to a plunge and a merger. Though some aspects 
of this signal need three-dimensional numerical relativity simulations to be re- 
liably computed, we shall argue here for the need of developing in parallel an 
analytical description of the motion of coalescing binary black holes down to 
the merger phase. We shall present below an essential ingredient of such an 
analytical description: the 'effective one body ' approach to the motion of binary 
black holes [3|3]. 



2 Motion of binary pulsars in general relativity 

The traditional (text book) approach to the problem of motion of N separate 
bodies in General Relativity (GR) consists of solving, by successive approxima- 
tions, Einstein's field equations (we use the signature — h H — h) 

together with their consequence 

V V T» V = Q. (2) 
To do so, one assumes some specific matter model, say a perfect fluid, 

= {e+p)u^u v +pg^ . (3) 
One expands (say in powers of Newton's constant) 

9»v A ) = V»v + + h$ + . . . , (4) 
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and uses the simplifications brought by the 'Post-Newtonian' approximation 
(doh^ — c~ 1 d t h fl „ -C dih/iv] v/c Cl, p< e). Then one integrates the 
local material equation of motion @ over the volume of each separate body, 
labelled say by a = 1, 2, . . . , JV. In so doing, one must define some 'center of 
mass' z l a of body a, as well as some (approximately conserved) 'mass' m a of 
body a, together with some corresponding 'spin vector' S l a and, possibly, higher 
multipole moments. 

An important feature of this traditional method is to use a unique coor- 
dinate chart x^ to describe the full iV-body system. For instance, the center 
of mass, shape and spin of each body a are all described within this common 
coordinate system x^. This use of a single chart has several inconvenient as- 
pects, even in the case of weakly self-gravitating bodies (as in the solar system 
case). Indeed, it means for instance that a body which is, say, spherically sym- 
metric in its own 'rest frame' X a will appear as deformed into some kind of 
ellipsoid in the common coordinate chart x^ . Moreover, it is not clear how 
to construct 'good definitions' of the center of mass, spin vector, and higher 
multipole moments of body a, when described in the common coordinate chart 
x^. In addition, as we are interested in the motion of strongly self-gravitating 
bodies, it is not a priori justified to use a simple expansion of the type ^ be- 
cause h/jtu ~ ^2Gm a /(c 2 \x — z a \) will not be uniformly small in the common 

a 

coordinate system x^. It will be small if one stays far away from each object a, 
but, as recalled above, it will become of order unity on the surface of a compact 
body. 

These two shortcomings of the traditional 'one-chart' approach to the rela- 
tivistic problem of motion can be cured by using a 'multi-chart' approach.The 
multi-chart approach describes the motion of N (possibly, but not necessarily, 
compact) bodies by using N + l separate coordinate systems: (i) one global coor- 
dinate chart x^ (fx = 0, 1, 2, 3) used to describe the spacetime outside N 'tubes', 
each containing one body, and (ii) N local coordinate charts X" (a = 0, 1, 2, 3; 
a = 1,2,... ,N) used to describe the spacetime in and around each body a. 
The multi-chart approach was first used to discuss the motion of black holes 
and other compact objects [4l El O [7l O O [101 [H] . Then it was also found to 
be very convenient for describing, with the high-accuracy required for dealing 
with modern technologies such as VLBI, systems of N weakly self-gravitating 
bodies, such as the solar system [T2l [13] . 

The essential idea of the multi-chart approach is to combine the information 
contained in several expansions. One uses both a global expansion of the type 
((4]) and several local expansions of the type 

G aP {X2) = G^(X2; m a ) + H$(XJ;m a) m b ) + ■ ■ ■ , (5) 

where G^0(X;m a ) denotes the (possibly strong-field) metric generated by an 
isolated body of mass m a (possibly with the additional effect of spin). 

The separate expansions (|4]) and ([5]) are then 'matched' in some overlapping 
domain of common validity of the type Gm a /c 2 < R a <C \x — z a \ <C-D^ 
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I x a — x b \ (with b ^ a), where one can relate the different coordinate systems by 
expansions of the form 

x» = z£(T a ) + ef (T ) XI + \ /£(T ) X l a Xi + • • • (6) 

The multi-chart approach becomes simplified if one considers compact bodies 
(of radius R a comparable to 2 Gm a /c 2 ). In this case, it was shown [S], by 
considering how the 'internal expansion' ([5]) propagates into the 'external' one 
^ via the matching ([S]), that, in General Relativity, the internal structure of 
each compact body was effaced to a very high degree, when seen in the external 
expansion ([J). For instance, for non spinning bodies, the internal structure of 
each body (notably the way it responds to an external tidal excitation) shows 
up in the external problem of motion only at the fifth post-Newtonian (5PN) 
approximation, i.e. in terms of order (w/c) 10 in the equations of motion. 

This 'effacement of internal structure ' indicates that it should be possible 
to simplify the rigorous multi-chart approach by skeletonizing each compact 
body by means of some delta-function source. Mathematically, the use of dis- 
tributional sources is delicate in a nonlinear theory such as GR. However, it 
was found that one can reproduce the results of the more rigorous matched- 
multi-chart approach by treating the divergent integrals generated by the use 
of delta-function sources by means of (complex) analytic continuation [S] . The 
most efficient method (especially to high PN orders) has been found to use 
analytic continuation in the dimension of space d [14). 

Finally, the most efficient way to derive the general relativistic equations of 
motion of N compact bodies consists of solving the equations derived from the 
action (where g = — Aet{g^ v )) 

S = J d - X Vff R ^ -J2 maC j \J-9l^( z a) dz * dz a > (7) 

formally using the standard weak- field expansion (H|) , but considering the space 
dimension d as an arbitrary complex number which is sent to its physical value 
d = 3 only at the end of the calculation. This 'skeletonized' effective action 
approach to the motion of compact bodies has been extended to other theories 
of gravity [3 [11] . Finite-size corrections can be taken into account by adding 
nonminimal worldlinc couplings to the effective action (|Tj) [HI [16] . 

Using this method^ one has derived the equations of motion of two compact 
bodies at the 2.5PN (v 5 /c 5 ) approximation level needed for describing binary 
pulsars [13 ESI E]: 

d 2 z i 

-TjT = A l a0 {z a - z b )+c- 2 A l a2 {z a - z b ,v a ,v b ) 

+ c ~ iAt a4( Z a - Z b ,V a ,V b ,S a ,S b ) 

+ c- 5 Al 5 (z a -z b ,v a -v b )+0(c- 6 ). (8) 

1 Or, more precisely, an essentially equivalent analytic continuation using the so-called 
'Riesz kernels'. 
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Here A l aQ = —Gmb(z l a — z l b )/\z a — Zb\ 3 denotes the Newtonian acceleration, A l a2 
its 1PN modification, A l a4 its 2PN modification (together with the spin-orbit 
effects), and A l a5 the 2.5PN contribution of order v 5 /c 5 . [See the references 
above; or the review [19j . for more references and the explicit expressions of A2, 
A4 and A5.] It was verified that the term A r ab has the effect of decreasing the 
mechanical energy of the system by an amount equal (on average) to the energy 
lost in the form of gravitational wave flux at infinity. Note, however, that here 
A l a5 was derived, in the near zone of the system, as a direct consequence of the 
general relativistic propagation of gravity, at the velocity c, between the two 
bodies. This highlights the fact that binary pulsar tests of the existence of A l a5 
are direct tests of the reality of gravitational radiation. 

The 2.5PN equations of motion ([8]) are accurate enough for interpreting 
(together with the corresponding 'timing formula' discussed next) current and 
foreseeable binary pulsar data. In Section [6] below we shall discuss recent im- 
provements (3PN and 3.5PN) in the knowledge of the equations of motion and 
their use (in a suitably resummed form) for describing the last orbits of coalesc- 
ing binary black holes. 



3 Timing of binary pulsars in general relativity 

In order to extract observational effects from the equations of motion ([8]) one 
needs to go through two steps: (i) to solve the equations of motion ((8]) so as to 
get the coordinate positions Z\ and 22 as explicit functions of the coordinate 
time t, and (ii) to relate the coordinate motion z a {t) to the pulsar observables, 
i.e. mainly to the times of arrival of electromagnetic pulses on Earth. 

The first step has been accomplished, in a form particularly useful for dis- 
cussing pulsar timing, in Ref. [20]. There (see also [21]) it was shown that, 
when considering the full (periodic and secular) effects of the A2 ~ v 2 /c 2 terms 
in Eq. ([5]), together with the secular effects of the A4 ~ w 4 /c 4 and A$ ~ v 5 /c 5 
terms, the relativistic two-body motion could be written in a very simple 'quasi- 
Keplerian' form (in polar coordinates), namely: 

', dt + a = u — e t sin u , (9) 



9 - O = (1 + k) 2arctan 



1 + Gfl \ 2 U 

tan — 



1-eeJ 2 



(10) 



R = r ab = a R (l - e R cosw) , (11) 
r a = \z a - zcm\ = Or(l - e r cosn) , (12) 
r b = \z b - z C m\ = <v(l - e r > cosu) . (13) 

Here n = 2ir/Pb denotes the orbital frequency, k = A8/2ir = (&>)/n = 
(uj) -Ph/27r the fractional periastron advance per orbit, u an auxiliary angle ('rel- 
ativistic eccentric anomaly'), et, eg, en, e r and e r > various 'relativistic eccentric- 
ities' and an,a r and a r > some 'relativistic semi-major axes'. See [20j for the 
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relations between these quantities, as well as their link to the relativistic energy 
and angular momentum E, J. A direct study [21] of the dynamical effect of the 
contribution A$ ~ v 5 /c 5 in the equations of motion |5]) has checked that it led 
to a secular increase of the orbital frequency n(t) ~ n(0)+n(t — to), and thereby 
to a quadratic term in the 'relativistic mean anomaly' I = J ndt + a appearing 
on the left-hand side (L.H.S.) of Eq. ©: 



As for the contribution A4 ~ v 4 /c 4 it induces several secular effects in the 
orbital motion: various 2PN contributions to the dimensionless periastron pa- 
rameter k (S4 k ~ v 4 /c 4 + spin-orbit effects), and secular variations in the incli- 
nation of the orbital plane (due to spin-orbit effects). 

The second step in relating © to pulsar observations has been accomplished 
through the derivation of a 'relativistic timing formula' 22, 23] . The 'timing 
formula' of a binary pulsar is a multi-parameter mathematical function relating 
the observed time of arrival (at the radio-telescope) of the center of the N th 
pulse to the integer N. It involves many different physical effects: (i) dispersion 
effects, (ii) travel time across the solar system, (iii) gravitational delay due to 
the Sun and the planets, (iv) time dilation effects between the time measured 
on the Earth and the solar-system-barycenter time, (v) variations in the travel 
time between the binary pulsar and the solar-system barycenter (due to relative 
accelerations, parallax and proper motion), (vi) time delays happening within 
the binary system. We shall focus here on the time delays which take place 
within the binary system. [See Refs. [2U [25] for the use of timing effects linked 
to parallax and proper motion.] 

For a proper derivation of the time delays occurring within the binary sys- 
tem we need to use the multi-chart approach mentionned above. In the 'rest 
frame' (X® = cT a , X l a ) attached to the pulsar a, the pulsar phenomenon can be 
modelled by the secularly changing rotation of a beam of radio waves: 



where $ a is the longitude around the spin axis. [Depending on the precise defi- 
nition of the rest-frame attached to the pulsar, the spin axis can either be fixed, 
or be slowly evolving, see e.g. [13].] One must then relate the initial direction 
(©a, *&a), an d proper time T ai of emission of the pulsar beam to the coordinate 
direction and coordinate time of the null geodesic representing the electromag- 
netic beam in the 'global' coordinates x M used to describe the dynamics of the 
binary system [NB: the explicit orbital motion ([5])- (fT3"]) refers to such global 
coordinates x° = ct, x 1 ]. This is done by using the link ^ in which z l a denotes 
the global coordinates of the 'center of mass' of the pulsar, T a the local (proper) 
time of the pulsar frame, and where, for instance 



I ~ cr + n Q (t ~to) + - n(t - t ) 2 . 



(14) 




(15) 




(16) 
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Using the link © (with expressions such as (HHJ) for the coefficients ef , . . .) 
one finds, among other results, that a radio beam emitted in the proper direction 
N l in the local frame appears to propagate, in the global frame, in the coordinate 
direction n % where 

n l = N l + -- N* —— + O { J . (17) 

This is the well known 'aberration effect', which will then contribute to the 
timing formula. 

One must also write the link between the pulsar 'proper time' T a and the 
coordinate time t = x° /c — z®/c used in the orbital motion (l9|)- (fT3|) . This reads 

-t?dT2=g MV (a>;)dz£dzZ (18) 

where the 'tilde' denotes the operation consisting (in the matching approach) in 
discarding in g M „ the 'self contributions' ~ (Gm a /R a ) n , while keeping the effect 
of the companion (~ Gmb/r a b, etc. . .). One checks that this is equivalent (in 
the dimensional-continuation approach) to taking x^ — z% for sufficiently small 
values of the real part of the dimension d. To lowest order this yields the link 

/ 2Gm b v 2 a y ^ f / Gm b 1 v> a \ 

\ C 2 r ab C 2 J J \ C 2 r ab 2 C 2 / 

which combines the special relativistic and general relativistic time dilation 
effects. Hence, following [23] we can refer to them as the 'Einstein time delay'. 

Then, one must compute the (global) time taken by a light beam emitted 
by the pulsar, at the proper time T a (linked to Emission by (HH)), in the initial 
global direction n % (see Eq. (fl7|)). to reach the barycenter of the solar system. 
This is done by writing that this light beam follows a null geodesic: in particular 

= ds 2 = g^(x x ) dx" dx v ~ - (l - c 2 dt 2 + fl + dx 2 (20) 

where U = Gm a /\x — z a \ + Grnb/\x — z b \ is the Newtonian potential within the 
binary system. This yields (with t e = Emission, t a = Wival) 

t a ~t e = t dtc- f a \dx\ + 1 t (^^- + \dx\ . (21) 

Jt e C JU C Jte \\X-Z a \ \X-Z b \j 

The first term on the last RHS of Eq. (|2"Tj) is the usual 'light crossing time' 
\ |-2barycentcr(ia) — z a {t e )\ between the pulsar and the solar barycenter. It con- 
tains the 'Roemer time delay' due to the fact that z a (t e ) moves on an orbit. The 
second term on the last RHS of Eq. (|2 1 1) is the 'Shapiro time delay' due to the 
propagation of the beam in a curved spacetime (only the Gmb piece linked to 
the companion is variable). For a discussion of the 0{v/c) fractional corrections 
to the Shapiro time delay see [26j and references therein. 
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When inserting the 'quasi-Keplerian' form ©~(H3]) of the relativistic motion 
in the 'Roemer' term in (j2~lj) , together with all other relativistic effects, one finds 
that the final expression for the relativistic timing formula can be significantly 
simplified by doing two mathematical transformations. One can redefine the 
'time eccentricity' e t appearing in the 'Kepler equation' @, and one can define 
a new 'eccentric anomaly' angle: u — > u ncw [we henceforth drop the superscript 
'new' on u}. After these changes, the binary-system part of the general relativis- 
tic timing formula [23] takes the form (we suppress the index a on the pulsar 
proper time T a ) 



tbarycenter - t = -D -1 [T + A R (T) + A E {T) + A-j(T) + A A (T)] 



(22) 



with 



Afl = x sinw[cos u — e(l + 8 r )\ + x[l — e 2 (l + Sg) 2 } 1 ^ 2 cos u> sin u , (23) 

A E = 7sin-u, (24) 

A5 = — 2rln{l — ecosu — s[sincj(cosw — e) + (1 — e 2 ) 1 / 2 coswsinw]}{25) 

Aa — A{s'm[u> + A e (u)] + esinw} + B{cos[ui + A e (u)] + ecosw} , (26) 

where x — xq + x(T — To) represents the projected light-crossing time (x — 
Opuisar sini/c), e — eo + e(T — To) a certain (relativistically-defined) 'timing 
eccentricity', A e {u) the function 

A e {u) = 2 arctan 



f 



1/2 



u 

tan- 



(27) 



u> = ujq + k A e (u) the 'argument of the periastron', and where the (relativisti- 
cally-defined) 'eccentric anomaly' u is the function of the 'pulsar proper time' 
T obtained by solving the Kepler equation 



e sin u = 2tt 



T-Tn 



P b 2 V P b 



(28) 



It is understood here that the pulsar proper time T corresponding to the N th 
pulse is related to the integer N by an equation of the form 



N = c + v v T + l - v v T 2 + i i> p T 3 



(29) 



From these formulas, one sees that 80 (and 8 r ) measure some relativistic distor- 
tion of the pulsar orbit, 7 the amplitude of the 'Einstein time delay@ Ae, and 
r and s the range and shape of the 'Shapiro time delayll A5. Note also that 

2 The post-Keplerian timing parameter 7, first introduced in |22| . has the dimension of 
time, and should not be confused with the dimensionless post-Newtonian Eddington parameter 
7 PPN probed by solar-system experiments (see below). 

3 The dimensionless parameter s is numerically equal to the sine of the inclination angle i 
of the orbital plane, but its real definition within the PPK formalism is the timing parameter 
which determines the 'shape' of the logarithmic time delay Ag(T). 
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the dimensionless PPK parameter k measures the non-uniform advance of the 
periastron. It is related to the often quoted secular rate of periastron advance 
<jj = (duu/dt) by the relation k = 6jP\ } /2it. It has been explicitly checked that 
binary-pulsar observational data do indeed require to model the relativistic pe- 
riastron advance by means of the non-uniform (and non-trivial) function of u 
multiplying k on the R.H.S. of Eq. (J27D [2?pl . Finally, we see from Eq. ([25]) that 
Pb represents the (periastron to periastron) orbital period at the fiducial epoch 
To, while the dimensionless parameter Pb represents the time derivative of Pi, 
(at T ). 

Schematically, the structure of the DD timing formula ([22]) is 

tbarycontcr ~ t = F [T N] {p K }; {p PK }; {q PK }} , (30) 

where ibarycenter denotes the solar-system barycentric (infinite frequency) ar- 
rival time of a pulse, T the pulsar emission proper time (corrected for aberra- 
tion), {p K } — {Pb, To, eo, ojo, xq} is the set of Keplerian parameters, {p PK = 
k, 7, Pb, r, s, 5g, e, i} the set of separately measurable post- Keplerian parameters, 
and {q PK } = {S r , A, B, D} the set of not separately measurable post- Keplerian 
parameters [27] . [The parameter D is a 'Doppler factor' which enters as an 
overall multiplicative factor D^ 1 on the right-hand side of Eq. (f22l).] 

A further simplification of the DD timing formula was found possible. In- 
deed, the fact that the parameters {q PK } = {S r , A, B, D} are not separately 
measurable means that they can be absorbed in changes of the other param- 
eters. The explicit formulas for doing that were given in [23] and [27 : they 
consist in redefining e,x,Pb,Sg and S r . At the end of the day, it suffices to 
consider a simplified timing formula where {S r , A, B, D} have been set to some 
given fiducial values, e.g. {0,0,0, 1}, and where one only fits for the remaining 
parameters {p K } and {p PK }. 

Finally, let us mention that it is possible to extend the general parametrized 
timing formula (|30|) by writing a similar parametrized formula describing the ef- 
fect of the pulsar orbital motion on the directional spectral luminosity [d(energy) 
/c?(time) d(frequency) d(solid angle)] received by an observer. As discussed in 
detail in [27] this introduces a new set of 'pulse-structure post-Keplerian pa- 
rameters'. 

4 Phenomenological approach to testing rela- 
tivistic gravity with binary pulsar data 

As said in the Introduction, binary pulsars contain strong gravity domains and 
should therefore allow one to test the strong-field aspects of relativistic gravity. 
The question we face is then the following: How can one use binary pulsar data 
to test strong-field (and radiative) gravity? 

4 Alas this function is theory-independent, so that the non-uniform aspect of the periastron 
advance cannot be used to yield discriminating tests of relativistic gravity theories. 
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Two different types of answers can be given to this question: a phenomeno- 
logical (or theory-independent) one, or various types of theory- dependent ap- 
proaches. In this Section we shall consider the phenomenological approach. 

The phenomenological approach to binary-pulsar tests of relativistic gravity 
is called the parametrized post-Keplerian formalism [281 127] . This approach 
is based on the fact that the mathematical form of the multi-parameter DD 
timing formula ([50)) was found to be applicable not only in General Relativity, 
but also in a wide class of alternative theories of gravity. Indeed, any theory 
in which gravity is mediated not only by a metric field but by a general 
combination of a metric field and of one or several scalar fields ip^ will induce 
relativistic timing effects in binary pulsars which can still be parametrized by the 
formulas (f2"2"| - (|29[) . Such general 'tensor-multi-scalar' theories of gravity contain 
arbitrary functions of the scalar fields. They have been studied in full generality 
in |29j . It was shown that, under certain conditions, such tensor-scalar gravity 
theories could lead, because of strong-field effects, to very different predictions 
from those of General Relativity in binary pulsar timing observations [50ll5Tl[T5] . 
However, the point which is important for this Section, is that even when such 
strong-field effects develop one can still use the universal DD timing formula 
(|5D|) to fit the observed pulsar times of arrival. 

The basic idea of the phenomenological, parametrized post-Keplerian (PPK) 
approach is then the following: By least-square fitting the observed sequence 
of pulsar arrival times t^ to the parametrized formula (|30|) (in which TV is 
defined by Eq. (|29p which introduces the further parameters u p ,v pi i) p ) one can 
phenomenologically extract from raw observational data the (best fit) values of 
all the parameters entering Eqs. (|29|) and (|30l . In particular, one so determines 
both the set of Keplerian parameters {p K } = {Pb, To, eo, loq, Xq}, and the set of 
post-Keplerian (PK) parameters {p PK } — {k, 7, Pb, r, s, 8$, e, x}. In extracting 
these values, we did not have to assume any theory of gravity. However, each 
specific theory of gravity will make specific predictions relating the PK param- 
eters to the Keplerian ones, and to the two (a priori unknown) masses m a and 
mjy of the pulsar and its companion. [For certain PK parameters one must also 
consider other variables related to the spin vectors of a and b.] In other words, 
the measurement (in addition of the Keplerian parameters) of each PK param- 
eter defines, for each given theory, a curve in the (m a ,mb) mass plane. For any 
given theory, the measurement of two PK parameters determines two curves 
and thereby generically determines the values of the two masses m a and mb (as 
the point of intersection of these two curves). Therefore, as soon as one mea- 
sures three PK parameters one obtains a test of the considered gravity theory. 
The test is passed only if the three curves meet at one point. More generally, 
the measurement of n PK timing parameters yields n — 2 independent tests of 
relativistic gravity. Any one of these tests, i.e. any simultaneous measurement 
of three PK parameters can either confirm or put in doubt any given theory of 
gravity. 

As General Relativity is our current most successful theory of gravity, it is 
clearly the prime target for these tests. We have seen above that the timing data 
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of each binary pulsar provides a maximum of 8 PK parameters: k, 7, Pb, r, s,5g,e 
and x. Here, we were talking about a normal 'single line' binary pulsar where, 
among the two compact objects a and b only one of the two, say a is observed 
as a pulsar. In this case, one binary system can provide up to 8 — 2 = 6 tests of 
GR. In practice, however, it has not yet been possible to measure the parameter 
5g (which measures a small relativistic deformation of the elliptical orbit), nor 
the secular parameters e and x. The original Hulse- Taylor system PSR 1913+16 
has allowed one to measure 3 PK parameters: k = (dj)Pb/2ir, 7 and Pb- The 
two parameters k and 7 involve (non radiative) strong-field effects, while, as 
explained above, the orbital period derivative Pb is a direct consequence of the 
term A$ ~ v 5 jc° in the binary-system equations of motion ([5]). The term A$ is 
itself directly linked to the retarded propagation, at the velocity of light, of the 
gravitational interaction between the two strongly self-gravitating bodies a and 
b. Therefore, any test involving Pb will be a mixed radiative strong-field test. 

Let us explain on this example what information one needs to implement a 
phenomenological test such as the (k— 7— Pb)igi3 + i6 one. First, we need to know 
the predictions made by the considered target theory for the PK parameters k, 7 
and Pb as functions of the two masses m a and nib- These predictions have been 
worked out, for General Relativity, in Refs. [321[5TJ[33]. Introducing the notation 
(where n = 2-K/Pb) 

M = m a + nib (31) 
X a = m a /M ; X b = m b /M ; X a + X b = 1 (32) 



(GMn 



1/3 



they read 



MM) = — ^ — J , (33) 

fc GR (m Q ,m b ) = Po' ( 34 ) 

7 GR (m a ,m b ) = - X b {l + X b ) P 2 , (35) 
n 

1Q9tt 1 4- — p 2 -I- — p 4 

P b GR (m a ,m fe ) = -™f ( f_ e2 | 7/ 9 2 6 X a X h 0b- 06) 

However, if we use the three predictions (j34)) ~ ([3"6"|) . together with the best 
current observed values of the PK parameters fc obs , 7° bs , p° bd [32] we shall find 



b 

that the three curves fc GR (m a , nib) = fc obs , 7 GR (m a , nib) = 7° bs , -P GR (m a , nib ) — 
pobs - m (j ria ^ mb ^ mass plane jail to meet at about the 13 a level! Should 
this put in doubt General Relativity? No, because Ref. [33] has shown that the 
time variation (notably due to galactic acceleration effects) of the Doppler fac- 
tor D entering Eq. (|22|) entailed an extra contribution to the 'observed' period 
derivative P fc obs . We need to subtract this non-GR contribution before drawing 
the corresponding curve: P GR (m a , nib) — P & obs — ps alactK \ Then one finds that 
the three curves do meet within one a. This yields a deep confirmation of Gen- 
eral Relativity, and a direct observational proof of the reality of gravitational 
radiation. 
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We said several times that this test is also a probe of the strong-field aspects 
of GR. How can one see this? A look at the GR predictions (j34]) - (f36|) does not 
exhibit explicit strong-field effects. Indeed, the derivation of Eqs. ([3"4"]) - ([3"6"]) used 
in a crucial way the 'effacement of internal structure' that occurs in the general 
relativistic dynamics of compact objects. This non trivial property is rather 
specific of GR and means that, in this theory, all the strong-field effects can be 
absorbed in the definition of the masses m a and m?,. One can, however, verify 
that strong-field effects do enter the observable PK parameters k, 7, P& etc. . . by 
considering how the theoretical predictions ([3~4"]) - ([3l)]) get modified in alternative 
theories of gravity. The presence of such strong-field effects in PK parameters 
was first pointed out in Ref. [7] (see also [31]) for the Jordan-Fierz-Brans-Dickc 
theory of gravity, and in Ref. [8] for Rosen's bi-metric theory of gravity. A 
detailed study of such strong-field deviations was then performed in [29l [30l [31] 
for general tensor- (multi-)scalar theories of gravity. In the following Section 
we shall exhibit how such strong-field effects enter the various post-Keplerian 
parameters. 

Continuing our historical review of phenomenological pulsar tests, let us 
come to the binary system which was the first one to provide several 'pure 
strong-field tests' of relativistic gravity, without mixing of radiative effects: 
PSR 1534+12. In this system, it was possible to measure the four (non ra- 
diative) PK parameters k, 7, r and s. [We see from Eq. (|25|) that r and s mea- 
sure, respectively, the range and the shape of the 'Shapiro time delay' Ag.] The 
measurement of the 4 PK parameters fe,7, r, s define 4 curves in the (m a ,mf,) 
mass plane, and thereby yield 2 strong-field tests of GR. It was found in [35] 
that GR passes these two tests. For instance, the ratio between the measured 
value s obs of the phenomenological parameteJl s and the value s GR [/c obs , 7° bs ] 
predicted by GR on the basis of the measurements of the two PK parameters k 
and 7 (which determine, via Eqs. ([34]) , ([35]) . the GR-predicted value of m a and 
m b ) was found to be s obs /s GR [fc obs , 7° bs ] = 1.004 ±0.007 [35]. The most recent 
data [36] yield s ob 7s GR [fc obs , 7 obs ] = 1.000 ± 0.007. We see that we have here a 
confirmation of the strong-field regime of GR at the 1% level. 

Another way to get phenomenological tests of the strong field aspects of 
gravity concerns the possibility of a violation of the strong equivalence principle. 
This is parametrized by phenomenologically assuming that the ratio between the 
gravitational and the inertial mass of the pulsar differs from unity (which is its 
value in GR): (TO gra v/mi ne rt)a = 1+A Q . Similarly to what happens in the Earth- 
Moon-Sun system [37j . the three-body system made of a binary pulsar and of the 
Galaxy exhibits a 'polarization' of the orbit which is proportional to A = A a — 
Aj,, and which can be constrained by considering certain quasi-circular neutron- 
star- white-dwarf binary systems [38] . See [39j for recently published improved 
limitsll on the phenomenological equivalence-principle violation parameter A. 

5 As already mentioned the dimensionless parameter s is numerically equal (in all theories) 
to the sine of the inclination angle i of the orbital plane, but it is better thought, in the PPK 
formalism, as a phenomenological timing parameter determining the 'shape' of the logarithmic 
time delay Ag(T). 

6 Note, however, that these limits, as well as those previously obtained in |40| . assume 
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The Parkes multibeam survey has recently discovered several new interesting 
'relativistic' binary pulsars, thereby giving a huge increase in the number of 
phenomenological tests of relativistic gravity. Among those new binary pulsar 
systems, two stand out as superb testing grounds for relativistic gravity: (i) 
PSR J1141-6545 gUBK], and (ii) the remarkable double binary pulsar PSR 
J0737-3039A and B [H HQ Eg>] . 

The PSR J1141— 6545 timing data have led to the measurement of 3 PK 
parameters: k, 7, and Pb [42] , As in PSR 1913+16 this yields one mixed 
radiative-strong-held tesfl 

The timing data of the millisecond binary pulsar PSR J0737— 3039A have 
led to the direct measurement of 5 PK parameters: k, 7, r, s and Pf, [4"4"l 1431 
46 . In addition, the 'double line' nature of this binary system (i.e. the fact 
that one observes both components, A and B, as radio pulsars) allows one to 
perform new phenomenological tests by using Keplerian parameters. Indeed, the 
simultaneous measurement of the Keplerian parameters x a and Xb representing 
the projected light crossing times of both pulsars (^4 and B) gives access to the 
combined Keplerian parameter 

r obs 

R ohs = -±- . (37) 

X a 

On the other hand, the general derivation of |23J (applicable to any Lorentz- 
invariant theory of gravity, and notably to any tensor-scalar theory) shows that 
the theoretical prediction for the the ratio R, considered as a function of the 
masses m a and nib, is 

^theory = ?± +0 (!!L\ (38) 

m b \ c 4 J 

The absence of any explicit strong-held-gravity effects in the theoretical predic- 
tion (f3"8"]) (to be contrasted, for instance, with the predictions for PK parameters 
in tensor-scalar gravity discussed in the next Section) is mainly due to the con- 
vention used in [23] and [27) for defining the masses m a and nib- These are 
always defined so that the Lagrangian for two non interacting compact objects 



that the (a priori pulsar-mass dependent) parameter A ~ A a is the same for all the analyzed 
pulsars. 

7 In addition, scintillation data have led to an estimate of the sine of the orbital inclination, 
sine |47| . As said above, sini numerically coincides with the PK parameter s measuring the 
'shape' of the Shapiro time delay. Therefore, one could use the scintillation measurements as 
an indirect determination of s, thereby obtaining two independent tests from PSR J1141— 6545 
data. A caveat, however, is that the extraction of sini from scintillation measurements rests 
on several simplifying assumptions whose validity is unclear. In fact, in the case of PSR 
J0737— 3039 the direct timing measurement of s disagrees with its estimate via scintillation 
data 1401 . It is therefore safer not to use scintillation estimates of sin i on the same footing as 
direct timing measurements of the PK parameter s. On the other hand, a safe way of obtaining 
an s-related gravity test consists in using the necessary mathematical fact that s = s'mi < 1. 
In GR the definition x a = a, a sini/c leads to sini = nx a /(Po ^b)- Therefore we can write the 
inequality nx a /{0o{M) Xt,) < 1 as a phenomenological test of GR. 



13 



reads Lo = — m a c 2 (l — i^/c 2 ) 1 / 2 . In other words, m a c 2 represents the total 

a 

energy of body a. This means that one has implicitly lumped in the definition 
of m a many strong-self-gravity effects. [For instance, in tensor-scalar gravity 
m a includes not only the usual Einsteinian gravitational binding energy due 
to the self-gravitational field g^ u (x), but also the extra binding energy linked 
to the scalar field f(x).] Anyway, what is important is that, when performing 
a phenomenological test from the measurement of a triplet of parameters, e.g. 
{fc,7,i?}, at least one parameter among them be a priori sensitive to strong- 
field effects. This is enough for guaranteeing that the crossing of the three 
curves k thcm * {m ai m b ) = k ohs , 7 thcor y(TO a , m b ) = 7° bs , i? thcory (m a , m 6 ) = R obs 
is really a probe of strong-field gravity. 



PSRB1913+16 », PSRB1534+12 




0.5 1 1.5 2 2.5 m A 0.5 1 1.5 2 2.5 



PSRJ1 141-6545 m B PSR J0737-3039 




0.5 1 1.5 2 2.5 m A 0.5 1 1.5 2 2.5 m A 



Figure 1: Phenomenological tests of General Relativity obtained from Keplerian 
and post-Keplerian timing parameters of four relativistic pulsars. Figure taken 
from [48]. 

In conclusion, the two recently discovered binary pulsars PSR J1141— 6545 
and PSR J0737— 3039 have more than doubled the number of phenomenological 
tests of (radiative and) strong-field gravity. Before their discovery, the 'canoni- 
cal' relativistic binary pulsars PSR 1913+16 and PSR 1534+12 had given us four 
such tests: one (k — j — Pb) test from PSR 1913+16 and three (k — 7 — r — s — Pjl) 



8 The timing measurement of P° bs in PSR 1534+12 is even more strongly affected by 
kinematic corrections (D terms) than in the PSR 1913+16 case. In absence of a precise, 
independent measurement of the distance to PSR 1534+12, the k — 7 — P5 test yields, at best, 
a ~ 15% test of GR. 
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tests from PSR 1534+12. The two new binary systems have given us fiv^ more 
phenomenological tests: one (fc — 7 — P&) (or two, k — 7 — Pb — s) tests from PSR 
J1141-6545 and four (k--y-r-s-P b -R) tests from PSR J0737- 303^3- As 
illustrated in Figure [H these nine phenomenological tests of strong-field (and 
radiative) gravity are all in beautiful agreement with General Relativity. 

In addition, let us recall that several quasi-circular wide binaries, made of 
a neutron star and a white dwarf, have led to high-precision phenomenological 
confirmations |39j (in strong-field conditions) of one of the deep predictions of 
General Relativity: the 'strong' equivalence principle, i.e. the fact that var- 
ious bodies fall with the same acceleration in an external gravitational field, 
independently of the strength of their self-gravity. 

Finally, let us mention that Ref. [27] has extended the philosophy of the 
phenomenological (parametrized post-Keplerian) analysis of timing data, to a 
similar phenomenological analysis of pulse- structure data. Ref. [27] showed that, 
in principle, one could extract up to 11 'post-Keplerian pulse-structure param- 
eters'. Together with the 8 post-Keplerian timing parameters of a (single-line) 
binary pulsar, this makes a total of 19 phenomenological PK parameters. As 
these parameters depend not only on the two masses m a , m\, but also on the two 
angles A, r\ determining the direction of the spin axis of the pulsar, the maximum 
number of tests one might hope to extract from one (single-line) binary pulsar 
is 19 — 4 = 15. However, the present accuracy with which one can model and 
measure the pulse structure of the known pulsars has not yet allowed one to 
measure any of these new pulse-structure parameters in a theory-independent 
and model-independent way. 

Nonetheless, it has been possible to confirm the reality (and order of mag- 
nitude) of the spin-orbit coupling in GR which was pointed out [49], [50] to be 
observable via a secular change of the intensity profile of a pulsar signal. Confir- 
mations of general relativistic spin-orbit effects in the evolution of pulsar profiles 
were obtained in several pulsars: PSR 1913+16 [EQIHS], PSR B1534+12 [S3] 
and PSR J1141— 6545 [54]. In this respect, let us mention that the spin-orbit 
interaction affects also several PK parameters, either by inducing a secular evo- 
lution in some of them (see [27] ) or by contributing to their value. For instance, 
the spin-orbit interaction contributes to the observed value of the periastron 
advance parameter k an amount which is significant for the pulsars (such as 
1913+16 and 0737—3039) where k is measured with high-accuracy. It was then 
pointed out [55j that this gives, in principle, and indirect way of measuring the 
moment of inertia of neutron stars (a useful quantity for probing the equation 
of state of nuclear matter [56],I57|). However, this can be done only if one mea- 
sures, besides fc, two other PK parameters with 10 -5 accuracy. A rather tall 
order which will be a challenge to meet. 

The phenomenological approach to pulsar tests has the advantage that it can 

9 Or even six, if we use the scintillation determination of s in PSR J1141— 6545. 
10 The companion pulsar 0737— 3039B being non recycled, and being visible only during a 
small part of its orbit, cannot be timed with sufficient accuracy to allow one to measure any 
of its post-Keplerian parameters. 
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confirm or invalidate a specific theory of gravity without making assumptions 
about other theories. Moreover, as General Relativity has no free parameters, 
any test of its predictions is a potentially lethal test. From this point of view, it is 
remarkable that GR has passed with flying colours all the pulsar tests if has been 
submitted to. [See, notably, Fig. [T]] As argued above, these tests have probed 
strong-field aspects of gravity which had not been probed by solar-system (or 
cosmological) tests. On the other hand, a disadvantage of the phenomenological 
tests is that they do not tell us in any precise way which strong-field structures, 
have been actually tested. For instance, let us imagine that one day one specific 
PPK test fails to be satisfied by GR, while the others are OK. This leaves us in a 
quandary: If we trust the problematic test, we must conclude that GR is wrong. 
However, the other tests say that GR is OK. This example shows that we would 
like to have some idea of what physical effects, linked to strong-field gravity, 
enter in each test, or even better in each PK parameter. The 'effacement of 
internal structure' which takes place in GR does not allow one to discuss this 
issue. This gives us a motivation for going beyond the phenomenological PPK 
approach by considering theory-dependent formalisms in which one embeds GR 
within a space of alternative gravity theories. 

5 Theory-space approach to testing relativistic 
gravity with binary pulsar data 

A complementary approach to testing gravity with binary pulsar data consists 
in embedding General Relativity within a multi-parameter space of alternative 
theories of gravity. In other words, we want to contrast the predictions of GR 
with the predictions of continuous families of alternative theories. In so doing we 
hope to learn more about which structures of GR are actually being probed in 
binary pulsar tests. This is a bit similar to the well-known psycho-physiological 
fact that the best way to appreciate a nuance of colour is to surround a given 
patch of colour by other patches with slightly different colours. This makes it 
much easier to detect subtle differences in colour. In the same way, we hope to 
learn about the probing power of pulsar tests by seeing how the phenomeno- 
logical tests summarized in Fig. [1] fail (or continue) to be satisfied when one 
continuously deform, away from GR, the gravity theory which is being tested. 

Let us first recall the various ways in which this theory-space approach has 
been used in the context of the solar-system tests of relativistic gravity. 

5.1 Theory-space approaches to solar-system tests of rel- 
ativistic gravity 

In the quasi-stationary weak-field context of the solar-system, this theory-space 
approach has been implemented in two different ways. First, the parametrized 
post-Newtonian (PPN) formalism [58l EH [60l [37l EU [62l HH [63] describes 
many 'directions' in which generic alternative theories of gravity might dif- 
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for in their weak-field predictions from GR. In its most general versions the 
PPN formalism contains 10 'post-Einstein' PPN parameters, 7 = 7 PPN — l ["1 . 
= /? PPN — 1, £, ai, a2, CI3, Ci, C2, C3i C4- Each one of these dimensionless quanti- 
ties parametrizes a certain class of slow-motion, weak-field gravitational effects 
which deviate from corresponding GR predictions. For instance, 7 parametrizes 
modifications both of the effect of a massive body (say, the Sun) on the light 
passing near it, and of the terms in the two-body gravitational Lagrangian which 
are proportional to (Gm a mb/r a b) • (v a — vi,) 2 jc 2 . 

A second way of implementing the theory-space philosophy consists in con- 
sidering some explicit, parameter-dependent family of alternative relativistic 
theories of gravity. For instance, the simplest tensor-scalar theory of gravity 
put forward by Jordan [64j . Fierz |65j and Brans and Dicke [66) has a unique 
free parameter, say a 2 , = (2u B d + 3) _1 . When a 2 , — > 0, this theory reduces 
to GR, so that a 2 , (or 1/Vbd) measures all the deviations from GR. When 
considering the weak-field limit of the Jordan-Fierz-Brans-Dicke (JFBD) the- 
ory, one finds that it can be described within the PPN formalism by choosing 
7 = -2 a§(l + aly 1 , = and £ = on = Q = 0. 

Having briefly recalled the two types of theory-space approaches used to 
discuss solar-system tests, let us now consider the case of binary-pulsar tests. 



5.2 Theory-space approaches to binary-pulsar tests of rel- 
ativistic gravity 

There exist generalizations of these two different theory-space approaches to the 
context of strong-field gravity and binary pulsar tests. First, the PPN formalism 
has been (partially) extended beyond the 'first post-Newtonian' (1PN) order 
deviations from GR (<~ v 2 / c 2 + Gm / c 2 r) to describe 2PN order deviations from 

GR (j? + §t) ^ E3- Remarkably, there appear only two new parameters 

at the 2PN leveQ e and £■ Also, by expanding in powers of the self-gravity 
parameters of body a and b the predictions for the PPK timing parameters 
in generic tensor-multi-scalar theories, one has shown that these predictions 
depended on several 'layers' of new dimensionless parameters |29j . Early among 
these parameters one finds, the 1PN parameters 0,j and then the basic 2PN 
parameters e and (, but one also finds further parameters 0s, (00'), 0" , ■ ■ ■ 
which would not enter usual 2PN effects. The two approaches that we have just 
mentionned can be viewed as generalizations of the PPN formalism. 

There exist also useful generalizations to the strong-field context of the idea 
of considering some explicit parameter-dependent family of alternative theo- 

llr The PPN parameter 7 PPN is usually denoted simply as 7. To distinguish it from the 
Einstein-time-delay PPK timing parameter 7 used above we add the superscript PPN. In 
addition, as the value of 7 PPN in GR is 1, we prefer to work with the parameter 7 = 7 PPN — 1 
which vanishes in GR, and therefore measures a 'deviation' from GR in a certain 'direction' 
in theory-space. Similarly with (3 = /3 PPN — 1. 

12 When restricting oneself to the general class of tensor-multi-scalar theories. At the 1PN 
level, this restriction would imply that only the 'directions' 7 and j3 are allowed. 
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ries of relativistic gravity. Early studies [7] [HI GH] focussed either on the one- 
parameter JFBD tensor-scalar theory, or on some theories which are not con- 
tinuously connected to GR, such as Rosen's bimetric theory of gravity. Though 
the JFBD theory exhibits a marked difference from GR in that it predicts the 
existence of dipole radiation, it has the disadvantage that the weak field, solar- 
system constraints on its unique parameter a 2 , are so strong that they drastically 
constrain (and essentially forbid) the presence of any non-radiative, strong-field 
deviations from GR. In view of this, it is useful to consider other 'mini-spaces' 
of alternative theories. 

A two-parameter mini-space of theories, that we shall denot^H here as 
TM/?', /?")> was introduced in [2S]. This two-parameter family of tensor-bi-scalar 
theories was constructed so as to have exactly the same first post-Newtonian 
limit as GR (i.e. 7 = f3 = ■ ■ ■ = 0), but to differ from GR in its predictions 
for the various observables that can be extracted from binary pulsar data. Let 
us give one example of this behaviour of the T2 (/?',/?") class of theories. For 
a general theory of gravity we expect to have violations of the strong equiva- 
lence principle in the sense that the ratio between the gravitational mass of a 
self-gravitating body to its inertial mass will admit an expansion of the type 

m grav 1 

-p^^l + A a = l^- m c a + r ]2 cl + ... (39) 

where c a = —2 9 q^q measures the 'gravitational compactness' (or fractional 
gravitational binding energy, c a ~ — 2 Ef rw /m a c 2 ) of body a. The numerical 
coefficient 771 of the contribution linear in c a is a combination of the first post- 
Newtonian order PPN parameters, namely rji = 4/3 — 7 [37] . The numerical 
coefficient 1J2 of the term quadratic in c a is a combination of the 1PN and 2PN 
parameters. When working in the context of the T2(/?',/3") theories, the 1PN 
parameters vanish exactly {J3 — = 7) and the coefficient of the quadratic term 
becomes simply proportional to the theory parameter /?' : rj2 = \Bj3', where 
B w 1.026. This example shows explicitly how binary pulsar data (here the data 
constraining the equivalence principle violation parameter A = A a — see 
above) can go beyond solar-system experiments in probing certain strong-self- 
gravity effects. Indeed, solar-system experiments are totally insensitive to 2PN 
parameters because of the smallness of c a ~ Gm a /c 2 R a and of the structure 
of 2PN effects [67] . By contrast, the 'compactness' of neutron stars is of order 
c a ~ 0.21 m a /M & ~ 0.3 [29 so that the pulsar limit |A| < 5.5 x 10~ 3 [H] yields, 
within the T 2 (f3',f3") framework, a significant limit on the dimensionless (2PN 
order) parameter f3' : \(3'\ < 0.12. 

Ref. [3T] introduced a new two-parameter mini-space of gravity theories, 
denoted here as Ti(a 0l (3 ), which, from the point of view of theoretical physics, 
has several advantages over the T2 (/?',/?") mini-space mentionned above. First, 
it is technically simpler in that it contains only one scalar field ip besides the 



13 We add here an index 2 to T as a reminder that this is a class of tensor- 6i-scalar theories, 
i.e. that they contain two independent scalar fields tpi,tfi2 besides a dynamical metric 
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metric g M „ (hence the index 1 on Ti(ao,/?o))- Second, it contains only positive- 
energy excitations (while one combination of the two scalar fields of T2(/3',/3") 
carried negative-energy waves). Third, it is the minimal way to parametrize 
the huge class of tensor-mono-scalar theories with a 'coupling function' a((p) 
satisfying some very general requirements (see below). 

Let us now motivate the use of tensor-scalar theories of gravity as alternatives 
to general relativity. 



5.3 Tensor-scalar theories of gravity 

Let us start by recalling (essentially from [3T]) why tensor- (mono)-scalar theories 
define a natural class of alternatives to GR. First, and foremost, the existence 
of scalar partners to the graviton is a simple theoretical possibility which has 
surfaced many times in the development of unified theories, from Kaluza-Klein 
to superstring theory. Second, they are general enough to describe many inter- 
esting deviations from GR (both in weak-field and in strong field conditions), 
but simple enough to allow one to work out their predictions in full detail. 

Let us therefore consider a general tensor-scalar action involving a metric 
g^v (with signature 'mostly plus'), a scalar field $, and some matter variables 
ip m (including gauge bosons): 

S = [ — g 1/2 \f@)R - Z($)r v d^ d v $ - U{$j\ + S m [1> m ;g MV ] . 

107T Cx* J C L J 

(40) 

For simplicity, we assume here that the weak equivalence principle is satisfied, 
i.e., that the matter variables ipm are all coupled to the same 'physical metric! 14 ! 
g^v. The general model (|40|) involves three arbitrary functions: a function F($>) 
coupling the scalar $ to the Ricci scalar of g^, R = i?(g Ml/ ), a function Z(&) 
rcnormalizing the kinetic term of and a potential function [/(<£>). As we 
have the freedom of arbitrary redefinitions of the scalar field, <f> — > $' = /(<&), 
only two functions among F, Z and U are independent. It is often convenient 
to rewrite (I40|) in a canonical form, obtained by redefining both $ and g M „ 
according to 

g; v = F(<S>)g^, (41) 

1/2 



tp = ± d$ 



3 F' 2 ($) 1 Z($) 

4 F 2 ($) + 2 F($) 



(42) 



14 Actually, most unified models suggest that there are violations of the weak equivalence 
principle. However, the study of general string-inspired tensor-scalar models |68) has found 
that the composition-dependent effects would be negligible in the gravitational physics of 
neutron stars that we consider here. The experimental limits on tests of the equivalence 
principle would, however, bring a strong additional constraint of order 10 — 5 a§ ~ Aa/a < 
10 -12 . As this constraint is strongly model-dependent, we will not use it in our exclusion 
plots below. One should, however, keep in mind that a limit on the scalar coupling strength 
o'q of order Oq < 10 -7 1681 T69 is likely to exist in many, physically-motivated, tensor-scalar 
models. 
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This yields 

s = y^q- J ^ aJ 2 [R* - 2<?r^^W - vfr)] + s m [^a 2 ^) g ; u ] , 

(43) 

where i?* = Rig^), where the potential 

%) = r 2 ($)[/($), (44) 

and where the conformal coupling function A(ip) is given by 

Mf) = F~ 1/2 (<&) , (45) 

with Q((p) obtained by inverting the integral (14"2"1> . 

The two arbitrary functions entering the canonical form (|43[) are: (i) the con- 
formal coupling function A((p), and (ii) the potential function V(ip). Note that 
the 'physical metric' g^ (the one measured by laboratory clocks and rods) is 
conformally related to the 'Einstein metric' <7*„, being given by g^ v = A 2 (ip) t/*„. 
The canonical representation is technically useful because it decouples the two 
irreducible propagating excitations: the spin-0 excitations are described by ip, 
while the pure spin-2 excitations are described by the Einstein metric <7*„ (with 
kinetic term the usual Einstein-Hilbert action oc Rig*^))- 

In many technical developments it is useful to work with the logarithmic 
coupling function a((p) such that: 

a(ip) = In A(ip) ; A(ip) = e a{lp) . (46) 

In the case of the general model (|4"0|) this logarithmic^ coupling function is 
given by 

afo) = ~ lnF($) , 

where must be obtained from (|42[). 

In the following, we shall assume that the potential V(ip) is a slowly varying 
function of ip which, in the domain of variation we shall explore, is roughly 
equivalent to a very small mass term V((p) ~ 2 m^ip — (fio) 2 with mz, of cosmo- 
logical order of magnitude vr? = O(Hq), or, at least, with a range = m^ 1 
much larger than the typical length scales that we shall consider (such as the 
size of the binary orbit, or the size of the Galaxy when considering violations of 
the strong equivalence principle). Under this assumption^! the potential func- 
tion V(ip) will only serve the role of fixing the value of ip far from the system 
(to <p(r = 00) = ipo), and its effect on the propagation of ip within the system 

5 As we shall mostly work with a(ip) below, we shall henceforth drop the adjective 'loga- 
rithmic'. 

16 Note, however, that, as was recently explored in [701 1711 172] , a sufficiently fast varying 
potential V(ip) can change the tensor-scalar phenomenology by endowing tp with a mass term 
mjl = ^ d 2 V/dip 2 which strongly depends on the local value of ip and, thereby can get large 
in sufficiently dense environments. 
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will be negligible. In the end, the tensor-scalar phenomenology that we shall 
explore only depends on one function: the coupling function a(tp). 

Let us consider some examples to see what kind of coupling functions might 
naturally arise. First, the simplest case is the Jordan-Fierz-Brans-Dicke action, 
which is of the general type (|4"0")) with 



F(*) 



(47) 
(48) 



where wbd is an arbitrary constant. Using Eqs. (|4"2")) , (|4"5|) above, one finds that 
— 2 ao ip = In $ and that the (logarithmic) coupling function is simply 



a(ip) = ao ip + const. 



(49) 



where ao = =f(2cjbd + 3) 1 / 2 , depending on the sign chosen in Eq. (|42]) . Inde- 
pendently of this sign, one has the link 



2 wbd 



(50) 



Note that 2wbd + 3 must be positive for the spin-0 excitations to have the 
correct (non ghost) sign. 

Let us now discuss the often considered case of a massive scalar field having 
a nonminimal coupling to curvature 



S = 



16ttG 

This is of the form (H0| with 



S m [i> 



F($) = 1 + £<I> 2 , Z($) = 1 , U($) = m|$ 2 



m J ■ 

(51) 
(52) 



The case £ = — | is usually referred to as that of 'conformal coupling'. With the 
variables (J5TJ> the theory is ghost-free only if 2 (1 + £$ 2 ) 2 (dip/d§) 2 = 1 + £(1 + 
6 £) $ 2 is everywhere positive. If we do not wish to restrict the initial values of 
we must have £(1 + 6 £) > 0. Introducing then the notation \ = v£(l + 6£), 
we get the following link between <E> and y>: 



2V2y> = | In 
+ V6 In 



l + 2x$ v / l + X 2 * 2 + X* 



1 - 2\/6£$ 



1 + £$ 2 



(53) 



For small values of $, this yields <p = + <D(§ 3 ). The potential and the 

coupling functions are given by 



l + £$ 2 



(54) 
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a(cp) =-±ln(l + & 2 ). (55) 

These functions have singularities when 1 +£$ 2 vanishes. If we do not wish 
to restrict the initial value of $ we must assume £ > (which then implies our 
previous assumption £(l + 6£) > 0). Then there is a one-to-one relation between 
$ and ip over the entire real line. Small values of $ correspond to small values 
of ip and to a coupling function 

a(<p) = -Zip 2 + 0(tf i ). (56) 

On the other hand, large values of |4>| correspond to large values of \ip\, and to 
a coupling function of the asymptotic form 

c 

a(ip) ~ -V2 - \tp\ + const. (57) 
X 

The potential V(ip) has a minimum at tp = 0, as well as other minima at 
tp — ► ±oo. If we assume, for instance, that m\ and the cosmological dynamics 
are such that the cosmological value of tp is currently attracted towards zero, 
the value of tp at large distances from the local gravitating systems we shall 
consider will be tp Q <C 1. 

As a final example of a possible tensor-scalar gravity theory, let us discuss 
the string-motivated dilaton-runaway scenario considered in [73). The starting 
action (a functional of g^v and $) was taken of the general form 

\ a' a' 
_(V$) 2 ] - - B F ($)F 2 - + • 

and it was assumed that all the functions have a regular asymptotic be- 

havior when $ — > +00 of the form = Ci + 0(e~*). Under this assumption 

the early cosmological evolution can push $ towards +00 (hence the name 'run- 
away dilaton'). In the canonical, 'Einstein frame' representation (|43p. one has, 
for large values of <&,<!> ~ c (p, where c is a numerical constant, and the coupling 
function to hadronic matter is given by 



a ^ cx A QCD (^) cx B- 1 ' 2 ^) exp[-8^ 2 63- 1 B F (<p)} 



where 63 is the one-loop rational coefficient entering the renormalization-group 
running of the gauge field coupling g 2 F . This finally yields a coupling function 
of the approximate form (for large values of ip) : 

a(ip) ~ k eT cip + const. , 

where the dimensionless constants A; and c are both expected to be of order unity. 
[The constant c must be positive, but the sign of k is not a priori restricted.] 
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Summarizing: the JFBD model yields a coupling function which is a linear 
function of <p, Eq. (|49p . a nonminimally coupled scalar yields a coupling function 
which interpolates between a quadratic function of ip, Eq. (|56p . and a linear one, 
Eq. (|57[) , and the dilaton-runaway scenario of Ref. [73] yields a coupling function 
of a decaying exponential type. 



5.4 The role of the coupling function a(ip); definition of the 
two-dimensional space of tensor-scalar gravity theo- 
ries Ti(a , A)) 

Let us now discuss how the coupling function a(tp) enters the observable pre- 
dictions of tensor-scalar gravity at the first post-Newtonian (1PN) level, i.e., 
in the weak-field conditions appropriate to solar-system tests. It was shown in 
previous work that, if one uses appropriate units in the asymptotic region far 
from the system, namely units such that the asymptotic value a(<po) of a(ip) 
vanished, all observable quantities at the 1PN level depend only on the values 
of the first two derivatives of the a(tp) at ip — tpo- More precisely, if one defines 

. . daiw) „, . da((p) d 2 a(<p) ,„„,. 
a = — ^ ; P(<p) = — ^ = — , 58 
dip dip o<p z 

and denotes by ao = a(po), (3q = /3(</?o) their asymptotic values, one finds 
(see, e.g., that the effective gravitational constant between two bodies (as 
measured by a Cavendish experiment) is given by 

G = G.(l + a§), (59) 

while, among the PPN parameters, only the two basic Eddington ones, 7 = 
7 PPN — 1, and (3 = /3 PPN — 1, do not vanish, and are given by 

7^7 PPN -1 = -2^, (60) 
3 - /5 PPN -1 = 1 aoPoa ° (61) 

The structure of the results (l60|) and (|6T|) can be transparently expressed by 
means of simple (Feynman-like) diagrams (see, e.g., [Zl])- Eqs. ([59]) and ((60]) 
correspond to diagrams where the interaction between two worldlines (repre- 
senting two massive bodies) is mediated by the sum of the exchange of one 
graviton and one scalar particle. The scalar couples to matter with strength 
~ ao \fG„. The exchange of a scalar excitation then leads to a term cx a^. On 
the other hand, Eq. (|6ip corresponds to a nonlinear interaction between three 
worldlines involving: (i) the 'generation' of a scalar excitation on a first world- 
line (factor ao), (ii) a nonlinear vertex on a second worldline associated to the 



17 In these units the Einstein metric g*„ and the physical metric <? M „ asymptotically coin- 
cide. 
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quadratic piece of a(ip) (a qna< i(ip) = \ 0o((f — fo) 2 ; so that one gets a factor /3 ), 
and (iii) the final 'absorption' of a scalar excitation on a third worldline (second 
factor ao). 

Eqs. ([60|) and (jBTj) can be summarized by saying that the first two coefficients 
in the Taylor expansion of the coupling function a(ip) around (p = (po (after 
setting a(tpo) = 0) 

a {v) = a o{f — <Po) + ^ 0a(<P — Vo) 2 H (62) 

suffice to determine the quasi-stationary, weak- field (1PN) predictions of any 
tensor-scalar theory. In other words, the solar-system tests only explore the 
'osculating approximation' (|6"2"|) (slope and local curvature) to the function a(<p). 
Note that GR corresponds to a vanishing coupling function a(<p) = (so that 
a o — 0o = ■ ■ ■ = 0), the JFBD model corresponds to keeping only the first term 
on the R.H.S. of (|62|) , while, for instance, the nonminimally coupled scalar field 
(with asymptotic value <po <C 1) does indeed lead to nonzero values for both ao 
and (3o, namely 

a ^-2£vo; /?o^-2£. (63) 
Finally the dilaton-runaway scenario considered above leads also to non zero 
values for both ao and /3o, namely 

a ~-kc e~ cva ; (3 ~+kc 2 e~ CVo , (64) 

for a largish value of <pa . Note that the dilaton-runaway model naturally predicts 
that ao 1, and that (3q is of the same order of magnitude as ao : A) — — ca o 
with c being (positive and) of order unity. The interesting outcome is that such 
a model is well approximated by the usual JFBD model (with /3q = 0). This 
shows that a JFBD-like theory could come out from a model which is initially 
quite different from the usual exact JFBD theory. 

As we shall discuss in detail below, solar-system tests constrain a§ and a^ |/3o| 
to be both small. This immediately implies that |ao| must be small, i.e., that 
the scalar field is linearly weakly coupled to matter. On the other hand, the 
quadratic coupling parameter @q is not directly constrained. Both its magnitude 
and its sign can be more or less arbitrary. Note that there are no a priori sign 
restrictions on 0o. The conformal factor A 2 (ip) = exp(2a(iy9)) entering Eq. (|43|) 
had to be positive, but this leads to no restrictions on the sign of a(<p) and of 
its various derivative^. For instance, in the nonminimally coupled scalar field 
case, it seemed more natural to require £ > 0, which leads to a negative /3q in 
view of Eq. (|63)) . 

Let us summarize the results above: (i) the most general tensor-scalar the- 
orj0is described by one arbitrary function a(np); and (ii) weak-field tests depend 

18 As explained above, we assume here the presence of a potential term V((p) to fix the 
asymptotic value ipo of ip. If the potential V(ip) is absent (or negligible), the 'attractor 
mechanism' of Refs. |75l (i* would attract ip to a minimum of the coupling function a(ip), 
thereby favoring a positive value of 

19 Under the assumption that the potential V(ip) is a slowly-varying function of ip, which 
modifies the propagation of ip only on very large scales. 
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only on the first two terms, parametrized by «o and A), in the Taylor expansion 
(|62p of a(tp) around its asymptotic value ipo . 

From this follows a rather natural way to define a simple mini space of tensor- 
scalar theories. It suffices to consider the two-dimensional space of theories, say 
Ti(ao,/?o), defined by the coupling function which is a quadratic polynomial in 
<p HUES], say 

a ao A,{f) = ao{ip - ip ) + - f3 (ip - tp ) 2 . (65) 

As indicated, this class of theories depends only on two parameters: «o and 
/3q. The asymptotic value tpo of (f does not count as a third parameter (when 
using the form (I65j) ) because one can always work with the shifted field (p = 
ip — (po, with asymptotic value <^o = and coupling function a ao ^ (Cp) = ao (p + 
h 0o (p 2 . Moreover, as already said, the asymptotic value a((po) of a(ip) has also 
no physical meaning, because one can always use units such that it vanishes (as 
done in (|65| ). 

Note also that an alternative way to represent the same class of theories is 
to use a coupling function of the very simple form 

a p(<P) = \P<P 2 > ( 66 ) 

but to keep the asymptotic value ipo as an independent parameter. This class 
of theories is clearly equivalent to Ti(oto, fio), Eq. (|65p . with the dictionary: 

a Q = /3<po, f3 = f3. 



5.5 Tensor-scalar gravity, strong-field effects, and binary- 
pulsar observables 

Having chosen some mini-space of gravity theories, we now wish to derive what 
predictions these theories make for the timing observables of binary pulsars. To 
do this we need to generalize the general relativistic treatment of the motion and 
timing of binary systems comprising strongly self-gravitating bodies summarized 
above. Let us recall that this treatment was based on a multi-chart method, 
using a matching between two separate problems: (i) the 'internal problem' 
considers each strongly self-gravitating body in a suitable approximately freely 
falling frame where the influence of its companion is small, and (ii) the 'external 
problem' where the two bodies are described as effective point masses which 
interact via the various fields they are coupled to. Let us first consider the 
internal problem, i.e., the description of a neutron star in an approximately 
freely falling frame where the influence of the companion is reduced to imposing 
some boundary conditions on the tensor and scalar fields with which it interacts 
[3 El [29l [301 [31] ■ The field equations of a general tensor-scalar theory, as derived 
from the canonical action (|43[) (neglecting the effect of V((f)) read 

r; u = 2 d u <p + svr g* (t; v - ~ T* g ;^j , (67) 

\3 g ,(p = -4ttG» a(<p)T* , (68) 
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where T* v = 2 c (g*)^ 1 ^ 2 SSm/Sg*^ denotes the material stress-energy tensor in 
'Einstein units', and a(ip) the (^-derivative of the coupling function, see Eq. ([58]) . 
All tensorial operations in Eqs. (|67[) and (f6"5)) are performed by using the Einstein 
metric g* v . 

Explicitly writing the field equations (|57|) and (f6"5| for a slowly rotating 
(stationary, axisymmetric) neutron star, labelled^! A, leads to a coupled set of 
ordinary differential equations constraining the radial dependence of <7*„ and 
ip [311 176] . Imposing the boundary conditions g*^ u — > T]^ v , (p — > (p a at large 
radial distances, finally determines the crucial 'form factors' (in Einstein units) 
describing the effective coupling between the neutron star A and the fields to 
which it is sensitive: total mass mA(f a ), total scalar charge uiA^a), and inertia 
moment /^(Va)- As indicated, these quantities arc functions of the asymptotic 
value tp a of ip felt by the considered neutron staiF^I. They satisfy the relation 
LUAisPa) = ~dmA(<fia)/d(p a - From them, one defines other quantities that play 
an important role in binary pulsar physics, notably 

a A {>Pa) = = — 7. , 69 

m A oip a 

(70) 

as well as 

M^)--^. (71) 



The quantity a a, Eq. (|69[) . plays a crucial role. It measures the effective coupling 
strength between the neutron star and the ambient scalar field. If we formally let 
the self-gravity of the neutron A tend toward zero (i.e., if we consider a weakly 
self-gravitating object), the function aA^a) becomes replaced by a(<p a ) where 
a(ip) = da(ip)/d<p is the coupling strength appearing in the R.H.S. of Eq. (|6"51) . 
Roughly speaking, we can think of ctA^a) as a (suitable defined) average value 
of the local coupling strength a(ip(r)) over the radial profile of the neutron star 
A. 

It was pointed out in Refs. [30l E] that the strong self-gravity of a neutron 
star can cause the effective coupling strength aA(fa) to become of order unity, 
even when its weak-field counterpart ao = a{ip a ) is extremely small (as is im- 
plied by solar-system tests that put strong constraints on the PPN combination 



20 We henceforth use the labels A and B for the (recycled) pulsar and its companion, instead 
of the labels a and b used above. We henceforth use the label a to denote the asymptotic 
value of some quantity (at large radial distances within the local frame, X l A or Xg , of the 
considered neutron star A or B). 

21 This ip a is a combination of the cosmological background value </3o and of the scalar 
influence of the companion of the considered neutron star. It varies with the orbital period and 
is determined as part of the 'external problem' discussed below. Note that, strictly speaking, 
the label a (for asymptotic) should be indexed by the label of the considered neutron star: i.e. 
one should use a label a a (and a locally asymptotic value <fia A ) when considering the neutron 
star A, and a label as (with a corresponding Lpa B ) when considering the neutron star B. 
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Figure 2: Dependence upon the baryonic mass via of the coupling parameter 
a a in the theory Ti(ao,/?o) with ao — —0.014, /3q = —6. Figure taken from 

7= — 2 0.0/(1+0%)). This is illustrated, in the minimal context of the Ti(an, /?o) 
class of theories, in Figured] 

Note that when the baryonic mass rfiA of the neutron star is smaller than 
the critical mass rh cr ~ 1.24 M@ the effective scalar coupling strength a a of 
the star is quite small (because it is proportional to its weak-field limit ao = 
a(ip a )). By contrast, when rfiA > rn cr , \cxa\ becomes of order unity, nearly 
independently of the externally imposed a$ = a a = a(tp a ). This interesting 
non-perturbative behaviour was related in [30, 31 to a mechanism of spontaneous 
scalarization, akin to the well-known mechanism of spontaneous magnetization 
of ferromagnets. See also [48 for a simple analytical description of the behaviour 
of OtA- 

Let us also mention in passing that, in the case where A is a black hole, the 
effective coupling strength a a actually vanishes |29j . This result is related to 
the impossibility of having (regular) 'scalar hair' on a black hole. 

We have sketched above the first part of the matching approach to the mo- 
tion and timing of strongly self-gravitating bodies: the 'internal problem'. It 
remains to describe the remaining 'external problem'. As already mentionncd 
(and emphasized, in the present context, by Eardley [71 lllj). the most efficient 
way to describe the external problem is, instead of matching in detail the exter- 
nal fields (g£„, f) to the fields generated by each body in its comoving frame, to 
'skeletonize' the bodies by point masses. Technically this means working with 
the effective action 

S = J ^ gl /2 [R* - 2 <?r d,f d v f\ 

J2 c I mA(f(z A ))(-g; v (zA) dz>X dz» A ) 1/2 , (72) 

A J 
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where the function rriA^p) hi the last term on the R.H.S. is the function mA(<Pa) 
obtained above by solving the internal problem. Eq. ([72]) indicates that the ar- 
gument of this function is taken to be ip a — ip(za), i-e., the value that the scalar 
field (as viewed in the external problem) takes at the location z^ A of the center of 
mass of body A. However, as body A is described, in the external problem, as a 
point mass this causes a technical difficulty: the externally determined field ip{x) 
becomes formally singular at the location of the point sources, so that <p{za) is a 
priori undefined. One can either deal with this problem by coming back to the 
physically well-defined matching approach (which shows that <p(za) should be 
replaced by ip a , the value of <p in an intermediate domain Ra -C r <C \za — %b \ ) , 
or use the efficient technique of dimensional regularization. This means that the 
spacetime dimension D in Eq. ([72]) is first taken to have a complex value such 
that <p(za) is finite, before being analytically continued to its physical value 
D = 4. 



One then derives from the action (|T2[) two important consequences for the 
motion and timing of binary pulsars. First, one derives the Lagrangian de- 
scribing the relativistic interaction between TV strongly self-gravitating bodies 
(including orbital - (v/c) 2 effects, and neglecting C(w 4 /c 4 ) ones) [TT1 [7] 125} 154] . 
It is the sum of one-body, two-body and three-body terms. 

The one-body action has the usual form of the sum (over the label A) of the 
kinetic term of each point mass: 

LT~ h ° dy = -m A c 2 ^l-v\/c 2 

= -m A c 2 + -m A v 2 A + -mA^- + f -jj . (73) 

Here, we use Einstein units, and the inertial mass ttia entering Eq. (|73[) is rriA = 
mA^o), where ifo is the asymptotic value of tp far away from the considered 
TV-body system. 

The two-body action is a sum over the pairs A, B of a term L^ ody which 
differs from the GR-predicted 2-body Lagrangian in two ways: (i) the usual 
gravitational constant G appearing as an overall factor in i^ ody must be re- 
placed by an effective (body-dependent) gravitational constant (in the appro- 
priate units mentioned above) given by 

Gab = G4l + a A a B ), (74) 

and (ii) the relativistic (0(v 2 / c 2 )) terms in L^ ody contain, in addition to those 
predicted by GR, new velocity-dependent terms of the form 

X1T 2-body /_ G A Bm, A m B [v A ~ v B ) 2 

ab -(iab) — -a , (75) 

rAB c 

with 

IAB = IAB - 1 = -2— . (76) 

1 + OLA OLB 
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In these expressions a a = ola{wo) = dlnmA((po)/d<Po (see Eq. (jB"9"|) with ip a — ► 

Finally, the 3-body action is a sum over the pairs B, C and over ^4 (with 
A ^ B, A ^ C, but the possibility of having B = C) of 

r 3-body _ /-, . 7,A sGABGACmAm B mc , , 

abc —-\ l + *Pbc) 9 ~ ~z y") 



where 



Cf TAB TAG 



3A _ r a « _ 1 «b /3a ac 

2 (1 + a A as)(l + aA ac) 
with = daA(<Po)/d<po (see Eq. ([TO]) with </? a — > ip ). 

When comparing the strong-field results ([73]). (|7S|) . fTSjl to their weak-field 
counterparts (f55|) , ([50)1 , ([CTjl one sees that the body-dependent quantity 
replaces the weak-field coupling strength ao in all quantities which are linked 
to a scalar effect generated by body A. Note also that, in keeping with the 
'3-body' nature of Eq. (|77|) . the quantity j3g C — 1 is linked to scalar interactions 
which are generated in bodies B and C and which nonlinearly interact on body 
A. The notation used above has been chosen to emphasize that jab and [3g C 
are strong-field analogs of the usual Eddington parameters 7 PPN , /3 PPN , so that 
jab and P B q are strong-field analogs of the 'post-Einstein' 1PN parameters j 
and (3 (which vanish in GR). Indeed the usual PPN results for the post-Einstein 
terms in the C(l/c 2 ) 2-body and 3-body Lagrangians are obtained by replacing 
in Eqs. (75]) and (77]) jab -> 7, Bbc _> P and G ab -> G. 

The non-perturbative strong-field effects discussed above show that the strong 
self-gravity of neutron stars can cause jab and /3g C to be significantly different 
from their GR values 7 GR = 1, /3 GR = 1, in some scalar-tensor theories having 
a small value of the basic coupling parameter ao (so that -y PPN — 1 oc and 
/3 PPN — 1 oc Po «q are both small). For instance, Fig. [2] shows that it is possible 
to have aA ~ ct B ~ ± 0.6 which implies jab — 1 ~ — 0.53, i.e., a 50% deviation 
from GR! Even larger effects can arise in /3g C — 1 because of the large values 
that Pa — dctA/difo can reach near the spontaneous scalarization transition 

CHI 

Those possible strong-field modifications of the effective Eddington param- 
eters jab, Pbc w hich parametrize the 'first post-Keplerian' (1PK) effects 
(i.e., the orbital effects ~ v 2 /c 2 smaller than those entailed by the Lagrangian 
Y]^ rriA v 2 A + | ^ Gab m A m s I tab), can then significantly modify the usual 

A A^B 

GR predictions relating the directly observable parametrized post-Keplerian 
(PPK) parameters to the values of the masses of the pulsar and its compan- 
ion. As worked out in Rcfs. [TTJ [371 131] one finds the following modified 
predictions for the PPK parameters k = (uj)/n, r and s: 



k th (mA,m B ) = 



Gab (mA + m B )n\ 2/Z 



1-e 2 

l-\a A a B X A (3 B a\ + X B Pa a% 



1 + a.AOtB 6 (1 + <X4 «b) 2 



(79) 
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r (mA,m B ) = G Bm B , (80) 



s th (m A ,m B ) 



Xr 



Gab{vtia + tub) n 



-1/3 



(81) 



Here, the label A refers to the object which is timed ('the pulsar! 22 !), the label B 
refers to its companion, xa = aAsirii/c denotes the projected semi- major axis 
of the orbit of A (in light seconds), Xa = rn A / {fnA + ms) and Xb = tub / '(wa + 
ma) = 1 — Xa the mass ratios, n = 2it / P& the orbital frequency and Gob = 
G*(l + aotts) the effective gravitational constant measuring the interaction 
between B and a test object (namely electromagnetic waves on their way from 
the pulsar toward the Earth). In addition one must replace the unknown bare 
Newtonian G* by its expression in terms of the one measured in Cavendish 
experiments, i.e., G* = G/(l + a 2 ,) as deduced from Eq. ([59]) . 

The modified theoretical prediction for the PPK parameter 7 entering the 
'Einstein time delay' Ab, Eq. is more complicated to derive because one 

must take into account the modulation of the proper spin period of the pulsar 
caused by the variation of its moment of inertia I a under the (scalar) influence 
of its companion [IT] [7] [31] . This leads to 

th / \ e X B ( ' G A B{mA + m B )n^ 1 " 

Y {mA,m B ) - 1 



n 1 + a a cub V c3 
[X B (l + a A a B ) + l + kAOtB], (82) 

where fe^t^o) = — 91n/^(</jg) / 'd'Po ( see Eq. ([TT]) with ip a — > ip ). Numerical 
studies |31j show that \za can take quite large values. Actually, the quantity 
\za olb entering (|82p blows up near the scalarization transition when ao — > 
(keeping /?o < fixed) . In other words a theory which is closer to GR in weak- 
field conditions predicts larger deviations in the strong-field regime. 

The structure dependence of the effective gravitational constant Gab , Eq. (f74|) , 
has also the consequence that the object A does not fall in the same way 
as B in the gravitational field of the Galaxy. As most of the mass of the 
Galaxy is made of non strongly-self-gravitating bodies, A will fall toward the 
Galaxy with an acceleration oc Gao, while B will fall with an acceleration 
oc Gbo- Here, as above, Gao — Goa — G*(l + a Q a.A) is the effective gravi- 
tational constant between A and any weakly self-gravitating body. As pointed 
out in Ref. [38] this possible violation of the universality of free fall of self- 
gravitating bodies can be constrained by using observational data on the class 
of small-eccentricity long-orbital-period binary pulsars. More precisely, the 
quantity which can be observationally constrained is not exactly the violation 
&-AB = (Goa — Gqb)/G = (1 + ag) _1 (ao ola — ao«B) of the strong equivalence 
principle [which simplifies to A^o = (Goa — G)/G = (1 + aQ)~ 1 (a aA — «o) 



22 In the double binary pulsar, both the first discovered pulsar and its companion are pulsars. 
However, the companion B is a non recycled, slow pulsar whose motion is well described by 
Keplerian parameters only. 
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in the case of observational relevance where one neglects the self-gravity of the 
white-dwarf companion ] but rathe@ [29] 

effective = (2 lAB -{X A ffi A +X B Pj SB )+2 



3 

(1 + a A a B r 3/2 (l + alTHao a A - a a B ) ■ (83) 

Here, the index B (= white-dwarf companion) can be replaced by (weakly self- 
gravitating body) so that, for instance, j AB = j A q — 1 — 2 a A ao/(l + ct A «o) = 
(1 — a A <*o)/(l + a A ot-o), as deduced from Eq. (|75|) . 

It remains to discuss the possible strong-field modifications of the theoretical 
prediction for the orbital period derivative Pb = P b h (m A ,m B ). This is obtained 
by deriving from the effective action (|72")l the energy lost by the binary system 
in the form of fluxes of spin-2 and spin-0 waves at infinity. The needed results 
in a generic tensor-scalar theory were derived in Refs. [29] [34] (in addition one 
must take into account the tensor-scalar modification of the additional 'varying- 
Doppler' contribution to the observed P b due to the Galactic acceleration [33] )• 
The final result for P b is of the form 

P'lth/ \ Amonopolc . Adipolc . Aquadrupolc . Aquadrupolc 

b (m A ,m B ) - r bip + r hlp + r blp + r bg , 

+ A S GR tiC + <* th p t laCtiC - (84) 



where, for instance, p™ ono P olc j s (heuristicallj{3) related to the monopolar flux 



nfi: 

quadrupolar flux of spin-2 waves at infinity. It reads: 



of spin-0 waves at infinity. The term pq^ adm P olc corresponds to the usual 



Aquadrupolo/ x 1927T m A m B 

P bg , p (m A ,m B ) = -77— 77 ^ (85) 

9 5(1 + a A a B ) (m A + m B y 

fG AB (m A + m B )n\ 5/3 1 + 73 e 2 /24 + 37 e 4 /96 



V c 3 J (1-e 2 ) 7 / 2 

with G AB = G*(l + a A a B ) = G(l + a A a B )/(l + 0$), where G» is the 'bare' 
gravitational constant appearing in the action, while G is the gravitational con- 
stant measured in Cavendish experiments. The flux (|85j) is the only one which 
survives in GR (although without any a^-related modifications). Among the 
several other contributions which arise in tensor-scalar theories, let us only write 
down the explicit expression of the contribution to (]84p coming from the dipolar 
flux of scalar waves. Indeed, this contribution is, in most cases, the dominant 



23 This refinement is given here for pedagogical completeness. However, in practice, the 
lowest-order result A ~ (1 + <^o) _ 1 ( a a A — <Xn) — ceo — oCq is accurate enough. 

24 Contrary to the GR case where a lot of effort was spent to show how the observed Pj, 
was directly related to the GR predictions for the (t)/c) 5 -accurate orbital equations of motion 
of a binary system [5], we use here the indirect and less rigorous argument that the energy 
flux at infinity should be balanced by a corresponding decrease of the mechanical energy of 
the binary system. 
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one [7J because it scales as (v/c) 3 , while the monopolar and quadrupolar con- 
tributions scale as {v/c) 5 . It reads 

ndipole, x G* m A m B n 1 + e 2 /2 2 

P bJ (mA,m B ) = -27T-T7 : c 7; 5^ "a - «s ■ (86) 



Note that the dipolar effect (|86|) vanishes when a a — as- Indeed, a binary 
system made of two identical objects (A — B) cannot select a preferred direction 
for a dipole vector, and cannot therefore emit any dipolar radiation. This also 
implies that double neutron star systems (which tend to have vtla ~ mj ~ 
1.35 Mq) will be rather poor emitters of dipolar radiation (though <(86j) still tends 
to dominate over the other terms in ([84| , because of the remaining difference 
(m,A — niB)/(mA + tub) 7^ 0). By contrast, very dissymmetric systems such 
as a neutron-star and a white-dwarf (or a neutron-star and a black hole) will 
be very efficient emitters of dipolar radiation, and will potentially lead to very 
strong constraints on tensor-scalar theories. See below. 

5.6 Theory-space analyses of binary pulsar data 

Having reviewed the theoretical results needed to discuss the predictions of 
alternative gravity theories, let us end by summarizing the results of various 
theory-space analyses of binary pulsar data. 

Let us first recall what are the best, current solar-system limits on the two 
1PN 'post-Einstein' parameters 7 = 7 PPN — 1 and (3 = /3 PPN — 1. They are: 

7= (2.1 ±2.3) x 1(T 5 , (87) 

from frequency shift measurements made with the Cassini spacecraft [75 , which 
supersedes the constraint 

7 = (-1.7 ±4.5) x 1CT 4 (88) 

from VLBI measurements [75] . 

|2 7 -/3| < 3 x 10~ 3 , (89) 

from Mercury's perihelion shift [531 HO] > and 

4/3-7 = (4-4 ± 4.5) x 1CT 4 , (90) 

from Lunar laser ranging measurements [81] . 

Concerning binary pulsar data, we can make use of the published measure- 
ments of various Keplerian and post-Keplerian timing parameters in the binary 
pulsars: PSR 1913+16 [3J, PSR B1534±12 [35], PSR J1141-6545 [Hj and 

PSR J0737-3039A±B [H [45l [46] . In addition . we can US S the recently up- 
dated limit on the parameter A measuring a possible violation of the strong 



25 There is, however, a caveat in the theoretical use one can make of the phenomenological 
limits on A. Indeed, in the small-eccentricity long-orbital-period binary pulsar systems used 
to constrain A one does not have access to enough PK parameters to measure the pulsar mass 
m,A directly. As the theoretical expression of A ~ ceo a A — c*q depends on itia (through a a), 
one needs to assume some fiducial value of ttia (say niA — 1-35 Mq). 
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equivalence principle (SEP), namely |A| < 5.5 x 10 3 at the 95% confidence 
level [39]. 

This ensemble of solar-system and binary-pulsar data can then be analyzed 
within any given parametrized theoretical framework. For instance, one might 
work within 

(i) the 4-parameter framework To (7, 0; e, £) [57] which defines the 2PN exten- 
sion of the original (Eddington) PPN framework To (7, 0); or 

(ii) the 2-parameter class of tensor-mono-scalar theories Ti(ao,/3o) [2Q]; or 

(iii) the 2-parameter class of tensor-bi-scalar theories T2(/3',/3") [2"9"] . 

Here, the index on To (7, (3; e, £) is a reminder of the fact that this framework 
is not a family of specific theories (it contains zero explicit dynamical fields), 
but is a parametrization of 2PN deviations from GR. As a consequence, its use 
for analyzing binary pulsar data is somewhat ill-defined because one needs to 
truncate the various timing observables (which are functions of the compactness 
of the two bodies A and B, say P PK = f(cA,cs)) at the 2PN order (i.e. es- 
sentially at the quadratic order in ca and/or eg). For some observables (or for 
product of observables) there might be several ways of defining this truncation. 
In spite of this slight inconvenience, the use of the To(j,(3;e,£) framework is 
conceptually useful because it shows very clearly why and how binary-pulsar 
data can probe the behaviour of gravitational theories beyond the usual 1PN 
regime probed by solar-system tests. 

For instance, the parameter = m^ rav /m™ crt — 1 measuring the strong 
equivalence principle (SEP) violation in a neutron star has, within the To (7, (3; e, Q 
framework, a 2PN-order expansion of the form [ST] 

A A = -i(4/3-7)c A +(!+C + 0(/3)) b A , (91) 

where c A = -2 ~ £ (U) A , b A = £ {U 2 ) A =2 B c\, with B ~ 1.026 and 

ca — kniA/MQ with k ~ 0.21. The general result (|91|) is compatible with the 
result quoted in subsection 15.21 within the context of the theory T 2 (/3',/3") when 
taking into account the fact that, within T 2 (/3',/3"), one has /3 = 7 = 0, e = fi' 
and C = [and that (3" parametrizes some effects beyond the 2PN level]. 

On the example of Eq. (|9"Tj) one sees that, after having used solar-system tests 
to constrain the first contribution on the RHS to a very small value, one can 
use binary-pulsar tests of the SEP to set a significant limit on the combination 
i € + £ of 2PN parameters. Other pulsar data then yield significant limits on 
other combinations of the two 2PN parameters e and £■ The final conclusion 
is that binary-pulsar data allow one to set significant limits (around or better 
than the 1% level) on the possible 2PN deviations from GR (in contrast to 
solar-system tests which are unable to yield any limit on e and £) [67] . For a 
recent update of the limits on e and £, which makes use of recent pulsar data 
see gg. 
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Let us now briefly discuss the use of mini-space of theories, such as Ti(ao, AO 
or T 2 (/3',/3"), for analyzing solar-system and binary-pulsar data. The basic 
methodology is to compute, for each given theory (e.g. for each given values of 
a and A) if one chooses to work in the Ti(ao, Ao) theory space) a goodness- 
of-fit statistics X 2 ( a o> A)) measuring the quality of the agreement between the 
experimental data and the considered theory. For instance, when considering 
the timing data of a particular pulsar, for which one has measured several PK 
parameters pi (i — 1, . . . ,n) with some standard deviations c°^ s , one defines, 
for this pulsar 

n 

X> ,AO = min V«T 2 (pf c ° ry (a , A>5 m A , m B ) - pf s f , (92) 

i=l 

where 'min' denotes the result of minimizing over the unknown masses tua, tub 
and where p* hcoiy (ao, Ai; in a, ms) denotes the theoretical prediction (within 
Ti(ao,/3o)) f° r the PK observable pi (given also the observed values of the 
Keplerian parameters). 

The goodness-of-fit quantity x 2 ( a o, AO wu l reach its minimum Xmin f° r some 
values, say a™ m , /So 1111 , of ao and A)- Then, one focusses, for each pulsar, on the 
level contours of the function 

Ax 2 (a ,/3o) = X 2 (ao,/3o)-Xmi„- (93) 

Each choice of level contour (e.g. A^ 2 = 1 or A^ 2 = 2.3) defines a certain 
region in theory space, which contains, with a certain corresponding 'confidence 
level', the 'correct' theory of gravity (if it belongs to the considered mini-space of 
theories). When combining together several independent data sets (e.g. solar- 
system data, and different pulsar data) we can define a total goodness-of-fit 
statistics x 2 o t( a 0; AoX by adding together the various individual x 2 ( a 0i Ao)- This 
leads to a corresponding combined contour A x 2 o t( a 0i AO- 
Let us end by briefly summarizing the results of the theory-space approach 
to relativistic gravity tests. For detailed discussions the reader should consult 
Refs. [2H [311 EE [77] , and especially the recent update [35] which uses the 
latest binary-pulsar data. 

Regarding the two-parameter class of tensor-bi-scalar theories T2(/3', (3") the 
recent analysis [IB] has shown that the A x 2 (/3\ A") corresponding to the double 
binary pulsar PSR J0737— 3039 was defining quite a small elliptical allowed 
region in the (A', A") plane. By contrast the other pulsar data define much wider 
allowed regions, while the strong equivalence principle tests define (in view of 
the theoretical result A ~ 1 + i Bj3'(c A — c 2 B )) a thin, but infinitely long, strip 
|A| < est. in the (A', 0") plane. This highlights the power of the double binary 
pulsar in probing certain specific strong-field deviations from GR. 

Contrary to the T 2 (A , (3") tensor-bi-scalar theories, which were constructed 
to have exactly the same first post-Newtonian limit as GF@ (so that solar- 
system tests put no constraints on A an d A")> the class of tensor- mono-scalar 

26 However, this could be achieved only at the cost of allowing some combination of the two 
scalar fields to carry a negative energy flux. 
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theories Ti(ao,/3o) is such that its parameters «o and 00 parametrize both the 
weak-field 1PN regime (see Eqs. (|60l) and (|6Tj) above) and the strong-field regime 
(which plays an important role in compact binaries) . This means that each class 
of solar-system data (see Eqs. (|57 |) — above) will define, via a corresponding 
goodness-of-fit statistics of the type, say 

_ 2 / o \ I ^Cassini \ — 2 /-theory / „ o \ -Cassini\2 

Xca SS ini( a o,Po) = (<7 7 ) (7 y (a ,^ )-7 ) 

a certain allowed regiorFI in the (arjiA)) plane. As a consequence, the analysis 
in the framework of the Ti(ao,/3o) space of theories allows one to compare 
and contrast the probing powers of solar-system tests versus binary-pulsar tests 
(while comparing also solar-system tests among themselves and binary-pulsar 
ones among themselves). The result of the recent analysis [15] is shown in 
Figure [3] 



matter 

• ♦ 

l«ol 




Figure 3: Solar-system and binary-pulsar constraints on the two-parameter fam- 
ily of tensor-mono-scalar theories T\(ao, /3o)- Figure taken from [48] . 

In Fig. [H the various solar-system constraints ([87]) -([90 ]) are concentrated 
around the horizontal po axis. In particular, the high-precision Cassini con- 
straint is the lower small grey strip. The various pulsar constraints are labelled 
by the name of the pulsar, except for the strong equivalence principle constraint 
which is labelled SEP. Note that General Relativity corresponds to the origin 
of the (ao,/3o) plane, and is compatible with all existing tests. 

The global constraint obtained by combining all the pulsar tests would, to a 
good accuracy, be obtained by intersecting the various pulsar-allowed regions. 

27 Actually, in the case of the Cassini data, as it is quite plausible that the positive value of 
the published central value ^Cassini _ +2.1 X is due to unsubtracted systematic effects, 

we use (jCassmi _ 2.3 x 10 -5 but 7 Cassml = 0. Otherwise, we would get unreasonably strong 
lcr limits on Qq because tensor-scalar theories predict that 7 must be negative, see Eqs. H60H 
and <6lt . 
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One can then see on Fig. [3] that it would be comparable to the pre-Cassini 
solar-system constraints and that its boundaries would be defined successively 
(starting from the left) by 1913+16, 1141-6545, 0737-3039, 1913+16 again 
and 1141—6545 again. 

A first conclusion is therefore that, at the quantitative level, binary-pulsar 
tests constrain tensor-scalar gravity theories as strongly as most solar-system 
tests (excluding the exceptionally accurate Cassini result which constrains a§ 
to be smaller than 1.15 x 10~ 5 , i.e. |ao| < 3.4 x 10~ 3 ). A second conclusion 
is obtained by comparing the behaviour of the solar-system exclusion plots and 
of the binary-pulsar ones around the negative /3q axis. One sees that binary- 
pulsar tests exclude a whole domain of the theory space (located on the left 
of [3q < —4) which is compatible with all solar-system experiments (even when 
including the very tight Cassini constraint). This remarkable qualitative feature 
of pulsar tests is a direct consequence of the existence of (non-perturbative) 
strong-field effects which start developing when the product — (3$ ca (with ca 
denoting, as above, the compactness of the pulsar) becomes of order unity. 

6 Motion and radiation of binary black holes: 
post-Newtonian-expanded results 

In Section[2]we mentioned that the 2.5PN accurate equations of motion (O were 
sufficiently accurate to interpret binary pulsar observations. By contrast, the 
forthcoming observations of gravitational wave signals from inspiralling binary 
black holes (and also inspiralling binary neutron stars, or mixed black-holc- 
neutron-star systems) has posed to theorists the double challenge of: (i) deriv- 
ing more accurate equations of motion, and (ii) deriving accurate expressions 
for the waveforms emitted by inspiralling, and even coalescing, compact bina- 
ries. Indeed, the premier targets for LIGO/VIRGO/GEO are the waveforms 
emitted during the late inspiral phase of compact binaries, as well as during the 
subsequent 'plunge' and 'merger' phases. During these phases the basic PN ex- 
pansion parameter 7 e = GM/c 2 D ~ (i>orbitai/c) 2 ceases to be numerically very 
small, and starts approaching values of order unity. It might then seem hopeless 
to tackle the motion and radiation of such close binary systems by means of a 
PN-expansion-type analytical approach. However, there are two reasons why it 
is meaningful, and probably very useful, to tackle the motion of close compact 
binaries by an analytical approach. 

The first reason is the need to describe with high accuracy the phasing of 
the gravitational waveforms emitted during the inspiral phase (i.e. before the 
plunge and merger). During the inspiral phase, the PN expansion parameter 
7e stays most of the time significantly below 1, though it increases to reach 
values of order ~ i at the end of the inspiral. It is then a priori reasonable 
to expect that the expansions of interesting physical quantities in powers of 
7e will converge sufficiently rapidly during most of the inspiral to allow one 
to deduce physically meaningful results from a PN expansion truncated at a 
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large enough order. This motivated the efforts of several groups for deriving 
equations of motion more accurate than the 2.5PN level mentioned above, and 
for deriving correspondingly accurate gravitational-wave generation formalisms. 
These efforts will be briefly reviewed in this Section. 

The second reason is that there might be ways of improving the convergence 
of PN expansions by using resummation methods. Two such methods have been 
particularly studied: one based on Pade approximants, and the other one on 
a novel approach to the dynamics of compact binaries called the 'effective one 
body' approach. 

Before discussing (in the next Section) these two resummation techniques, 
let us briefly recall the state of the art in analytical approaches to the motion 
and radiation of binary black holej^l. 

Two different gravitational-wave generation formalisms have been developed 
up to a high PN accuracy: (i) the Blanchet-Damour-Iyer formalism [82l [83l [HH 
[85l [86l [87l [88] combines a multipolar post-Minkowskian (MPM) expansion in 
the exterior zone with a post-Newtonian expansion in the near zone; while 
(ii) the Will-Wiseman-Pati formalism [53 [901 GUI 122] uses a direct integration 
of the relaxed Einstein equations. These formalisms were used to compute 
increasingly accurate estimates of the gravitational waveforms emitted by in- 
spiralling binaries. These estimates include both normal, near-zone generated 
post-Newtonian effects (at the 1PN 83., 2PN [33 |gl |55] , and 3PN [MUM] 
levels), and more subtle, wave-zone generated (linear and non- linear) 'tail ef- 
fects' [SHI H21 ISFJ [EE]- However, technical problems arose at the 3PN level. 
The representation of black holes by 'delta-function' sources causes the appear- 
ance of dangerously divergent integrals in the 3PN multipole moments. The 
use of Hadamard (partie finie) regularization did not allow one to unambigu- 
ously compute the needed 3PN-accurate quadrupole moment. Only the use of 
the (formally) diffeomorphism-invariant dimensional regularization method (i.e. 
analytic continuation in the dimension of space d) allowed one to complete the 
3PN-level gravitational-radiation formalism [99] . 

In parallel with the development of 3PN-accurate gravitational radiation 
formalisms, several groups (notably Jaranowski-Schafer and Blanchet-Faye) ex- 
tended the PN-type computation of the equations of motion of binary black 
holes beyond the 2.5PN level recalled in Section [2] above. Here also, the rep- 
resentation of black holes by delta- function sources and the use of the (non 
diffeomorphism invariant) Hadamard regularization method led to ambiguities 
in the computation of the badly divergent integrals that enter the 3PN equations 
of motion |10Q[ 1101] , By contrast, the use of the (diffeomorphism invariant) di- 
mensional regularization method allowed one to complete the determination of 
the 3PN-level equations of motion [1021 1103] . They have also been derived by 



28 For simplicity, we shall phrase the results in the context of binary black holes. Actually, 
we use the 'effacement of internal structure' mentioned above to skeletonize the black holes by 
means of delta-function sources. This means that, apart from quadrupole-deformation effects, 
the results are also valid for binary systems comprising neutron stars. 
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an Einstein-Infeld-Hoffmann-type surface-integral approach [104) . The 3.5PN 
terms in the equations of motion are also known [92[ 11051 1106] . 

The works mentioned in this Section (see [107] for a detailed account and 
more references) finally lead to PN-expanded results for the motion and radia- 
tion of binary black holes. For instance, the equations of motion are given in a 
form which generalize those of Section [5] above, namely (a = 1, 2; i = 1, 2, 3) 



~ z ± = Al cons + A lRR , (94) 



dt 2 



where 

A cons = A + c- 2 A 2 + cr 4 A A + c- 6 A e , (95) 
denotes the 'conservative' 3PN-accurate terms, while 

a rr = c -5 A5 + c -i Ai ^ ( Q6 ) 

denotes the time-asymmetric contibutions, linked to 'radiation reaction'. 

On the other hand, if we consider for simplicity the inspiralling motion of 
a quasi-circular binary system, the essential quantity describing the emitted 
gravitational waveform is the phase <j) of the quadrupolar gravitational wave 
amplitude h(t) ~ a(t) cos(<f>(t) + S). PN theory allows one to derive several 
different functional expressions for the gravitational wave phase <p, as a function 
either of time or of the instantaneous frequency. For instance, as a function 
of time, 4> admits the following explicit expansion in powers of 9 = vc 3 (t c — 
t)/5GM (where t c denotes a formal 'time of coalescence', M = mi + m-i and 
v = mi m2 /M 2 ) 

ct>{t) =<p c - v- 1 9 5 / 8 (l + J2(a n + a' n In 9) , (97) 

with some numerical coefficients a n , a' n which depend only on the dimensionless 
(symmetric) mass ratio v = mi rri2/M 2 . The derivation of the 3.5PN-accurate 
expansion (|97)l uses both the 3PN-accurate conservative acceleration (|95|) and 
a 3.5PN extension of the (fractionally) lPN-accurate radiation reaction accel- 
eration (|96[) obtained by assuming a balance between the energy of the binary 
system and the gravitational- wave energy flux at infinity (see, e.g., |107j ). 



7 Motion and radiation of binary black holes: 
the Effective One Body approach 

The PN-expanded results briefly reviewed in the previous Section are expected 
to yield accurate descriptions of the motion and radiation of binary black holes 
during their inspiralling stage, say up to the moment where the PN expansion 
parameter 7 e = GM/c 2 D reaches the value ~ h where the orbital motion is 
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expected to become dynamically unstable ('last stable (circular) orbit' and be- 
ginning of a 'plunge' leading to the merger of the two black holes). One possible 
strategy for having a complete description of the motion and radiation of binary 
black holes, covering all the stages (inspiral, plunge, merger, ring-down), would 
then be to try to 'stitch together' PN-expanded analytical results describing the 
inspiral phase with 3d numerical results describing the end of the inspiral, the 
plunge, the merger and the ring-down of the final black hole. 

However, we wish to argue that it might be possible to do a better use of all 
the analytical information contained in the PN-expanded results (l94|) - (l96l) . The 
basic claim (first made in [2 [3]) is that the use of suitable resummation meth- 
ods should allow one to describe, by analytical toold^l. a sufficiently accurate 
approximation of the entire waveform, from inspiral to ring-down, including the 
non-perturbative plunge and merger phases. To reach such a goal, one needs to 
make use of several tools: (i) resummation methods, (ii) exploitation of the flex- 
ibility of analytical approaches, (iii) extraction of the non-perturbative informa- 
tion contained in various numerical simulations, (iv) qualitative understanding 
of the basic physical features which determine the waveform. 

Before coming to grasp with some of these issues, let us emphasize some 
conceptual aspects of this programme. Recently, an important breakthrough in 
numerical relativity [TUHl EM 11101 EH Q3S EH] has led to the computation 
of the gravitational waveform emitted during the late inspiral, plunge, merger 
and ring-down of equal-mass {y — j), non-spinning binary black holes. Some 
sample numerical simulations have also begun to explore the multi-parameter 

space of coalescing unequal-mass (0 < v < |), spinning (Si, Sy black hole 
binaries (see, e.g., [1 141 II 15] V In spite of the high computer power used in these 
simulations, the calculation of one waveform, corresponding to specific values of 
the continuous parameters (y, Si, 6\, tpi, 62, ^2) parametrizing the considered 
initial binary state, takes a long time. It would be therefore extremely useful, 
for detection purposes, to have in hand a (quasi-)analytical approach which 
would combine the crucial non-perturbative information that we can get from 
numerical simulations, with the rich perturbative information that has been 
acquired in many years of work on the theory of the motion and radiation 
of binary black holes. The claim here is that the Effective One Body (EOB) 
approach offers enough flexibility in its definition and implementation to be able 
to smoothly combine these two types of information. First results towards this 
goal are given in Refs. [1161 033 031 031 EM EH] ■ 

Let us start by discussing the first tool used in the EOB approach: the 
systematic use of resummation methods. Two such methods have been employed 
(and combined), and some evidence has been given that they do significantly 



29 Here we use the adjective 'analytical' for methods that solve explicit (analytically given) 
ordinary differential equations (ODE), even if one uses standard (Runge-Kutta-type) numeri- 
cal tools to solve them. The important point is that, contrary to 3d numerical relativity simu- 
lations, numerically solving ODE's is extremely fast, and can therefore be done (possibly even 
in real time) for a dense sample of theoretical parameters, such as orbital {y = m\ va^/M, . . .) 
or spin (a\ = Si/ Gra\ , Si , ipx , . . .) parameters. 
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improve the convergence properties of PN expansions. The first method is the 
use of Pade approximants. It has been shown in Ref. 122J that near-diagonal 
Pade approximants of the radiation reaction force jffi seemed to provide a good 
representation of T down to the last stable orbit (which is expected to occur 
when D ~ 6GM/c 2 , i.e. when j e ~ i). The second method is a novel approach 
to the dynamics of compact binaries, which constitutes the core of the Effective 
One Body (EOB) method. The EOB method was introduced in [2] [3], and was 
further extended to the 3PN level in [123] , and by including spin effects in |124j . 

For simplicity of exposition, let us first explain the EOB method at the 
2PN level. The starting point of the method is the 2PN-accurate Hamilto- 
nian describing (in Arnowitt-Deser-Misner-type coordinates) the conservative, 
or time symmetric, part of the equations of motion ([5]) (i.e. the truncation 
A cons = A + c- 2 A 2 + c-*Aa of Eq. say H 2PN ( gi - ga.P1.P2). By going 

to the center of mass of the system (pi + p 2 = 0), one obtains a PN-expanded 
Hamiltonian describing the relative motion, q = q\ — q 2 , p = pi = —p2- 

#2PN iVC (<Z,P) = H (q,p) + ^H 2 (q,p) + ±H i (q,p) , (98) 

where H (q,p) = j^P 2 + (with M = mi + m 2 and fi = mim 2 /M) 

corresponds to the Newtonian approximation to the relative motion, while H 2 
describes 1PN corrections and H4 2PN ones. It is well known that, at the New- 
tonian approximation, H (q,p) can be thought of as describing a 'test particle' 
of mass [i orbiting around an 'external mass' GM. The EOB approach is a 
general relativistic generalization of this fact. It consists in looking for an 'ex- 
ternal spacetime geometry' gf^(x x ; GM) such that the geodesic dynamics of a 
'test particle' of mass fi within g e ^{x x , GM) is equivalen%^\ (when expanded in 
powers of 1/c 2 ) to the original, relative PN-expanded dynamics (|98|) . The ad- 
vantage of the EOB method is that it corapactifies the information contained in 
the rather complicated PN-expanded Hamiltonian (f9"8"| into the much simpler 
PN-expansions of the two independent metric coefficients A(R), B(R) of the 
'external' geometry 

effi dx" dx u = —A(R) c 2 dT 2 + B(R) dR 2 + R 2 (d6 2 + sin 2 9 dip 2 ) . (99) 

For instance, the crucial 'ggo*' metric coefficient A(R) (which fully encodes the 
energetics of circular orbits) is originally given, at 2PN order, by the PN expan- 
sion 

A 2PN {R) = I -2u + 2isu 3 , (100) 

30 We henceforth denote by T the Hamiltonian version of the radiation reaction term A RR , 
Eq. l |96|l . in the (PN-expanded) equations of motion. It can be heuristically computed up to 
(absolute) 5.5PN 95, 125, 99 and even 6PN 130 order by assuming that the energy radiated 
in gravitational waves at infinity is balanced by a loss of the dynamical energy of the binary 
system. 

31 See the above references to see the precise sense in which the two dynamics are equivalent. 
Let us just say here that the best way to think about it is to think of both dynamics in quantum 
terms as two sets of quantized energy spectra E^^ lve {n,€), E ex t(n,£) that are required to 
be mapped onto each other by an energy rescaling i?2pflj tlvc = /(-Eoxt)- 
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where u = GM/c 2 R and v = fi/M = mi mil(rn\ + m.2) 2 . 

The dimensionless parameter v = fJ,/M varies between (in the test mass 
limit mi <C m 2 ) and i (in the equal-mass case mi = m 2 ). When v — > 0, 
Eq. (|100p yields back, as expected, the well-known Schwarzschild time-time 
metric coefficient — <7oo hw = 1 — 2u = 1 — 2GM/c 2 R. One therefore sees in 
Eq. pOOp the role of v as a deformation parameter connecting a well-known 
test-mass result to a non trivial and new 2PN result. It is also to be noted 
that the 1PN EOB result Aipn(-R) = 1 — 2u happens to be ^-independent, and 
therefore identical to ^4 Schw = 1 — 2u. This is remarkable in view of the many 
non-trivial ^-dependent terms in the 1PN relative dynamics. The physically real 
1PN ^-dependence happens to be fully encoded in the function f(E) mapping 
the two energy spectra which was found to be always given by a very simple 
result: 

E EO b = (l + \ §f ) . (101) 

Let us emphasize the remarkable simplicity of the 2PN result (|100[) . The 
2PN Hamiltonian (j9"8|) contains eleven rather complicated ^-dependent terms. 
After transformation to the EOB format, the dynamical information contained 
in these eleven coefficients gets compactified into the very simple additional con- 
tribution + 2vu 3 in A(i?), together with an equally simple contribution in the 
radial metric coefficient: (A (R) B(R)) 2pn = 1 — Qvu 2 . This compactification 
process is even more drastic when one goes to the next (conservative) post- 
Newtonian order: the 3PN level, i.e. additional terms of order C(l/c 6 ) in the 
Hamiltonian (|9"5|) . As mentioned above, the complete obtention of the 3PN 
dynamics has represented quite a theoretical challenge and the final, resulting 
Hamiltonian is quite complicated. Even after going to the center of mass frame, 
the 3PN additional contribution \ He(q,p) to Eq. (|9"8]) introduces eleven new 
complicated ^-dependent coefficients. After transformation to the EOB format 
[2 , these eleven new coefficients get 'compactified' into only three additional 
terms: (i) an additional contribution to A(R), (ii) an additional contribution 
to B(R), and (iii) a 0(p 4 ) modification of the 'external' geodesic Hamiltonian. 
For instance, the crucial 3PN ^qq* metric coefficient becomes 

-43Pn(^) = l-2u + 2^u 3 + a 4 i/u 4 , (102) 

where 

94 41 

a 4 = -r- - 7T 2 ~ 18.6879027 . (103) 

The fact that the 3PN coefficient 04 in the crucial 'effective radial potential' 
^3Pn(-R), Eq. (|102[) . is rather large and positive indicates that the ^-dependent 
nonlinear gravitational effects lead, for comparable masses {v ~ |), to a last 
stable (circular) orbit (LSO) which has a higher frequency and a larger binding 
energy than what a naive scaling from the test-particle limit (y — > 0) would 
suggest. Actually, the PN-expanded form (|102[) of A 3 pn(-R) does not seem to 
be a good representation of the (unknown) exact function Aeob(-R) when the 
(Schwarzschild-like) relative coordinate R becomes smaller than about 6GM/c 2 
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(which is the radius of the LSO in the test-mass limit). It was therefore sug- 
gested [123] to further resuixF^l A 3 pn (R) by replacing it by a suitable Pade (P) 
approximant. For instance, the replacement of ^3pn(-R) by 

Al(R) = Pi[A 3PN (R)} = 1 -^r^-j 3 ( 104 ) 

ensures that the v = \ case is smoothly connected with the v = limit. 

The use of (I104p was suggested before one had any (reliable) non-perturbative 
information on the binding of close black hole binaries. Later, a compari- 
son with some 'waveless' numerical simulations of circular black hole binaries 
[116] has given some evidence that (| 104|) is physically adequate. There it was 
also emphasized that, in principle, the comparison between numerical data and 
EOB-based predictions should allow one to determine the effect of the unknown 
higher PN contributions to Eq. (|102p . For intance, one can add a 4PN-like 
term + a$vu 5 in Eq. (| 1 02|) , and then Pade the resulting radial function, say 
A\ = -P4PI3PN + &5^w 5 ]. Comparing the predictions of A^fas] to numerical 
data might then determine what is the physically preferred 'effective' value of 
the unknown coefficient 0,5. This is an example of the useful flexibility of analyt- 
ical approaches: the fact that one can tap numerically-based, non-perturbative 
information to improve the EOB approach. 

As recently emphasized |120) . it is quite useful to tap the information con- 
tained in the Regge-Wheeler-Zerilli-type signals emitted by test-particles orbit- 
ing black holes (small v limit). The numerical methods needed to compute 
these non-perturbative phenomena [1261 1127} 1119] are much simpler than the 
ones needed in the comparable- mass case [y ~ j), but their results contain a 
lot of useful physical insights, which are relatively easy to explore. 

Let us finally sketch the overall structure of the EOB approach to the motion 
and radiation of binary black holes. Two of the basic elements are the EOB 
Hamiltonian Heob{q,P, Si, S2) and the radiation reaction force T{q,p, Si, S2). 
We have indicated here that the EOB approach has been generalized to the case 
of arbitrarily spinning black holes EES Q3H [128] . This leads to ODE's for the 
evolution of the variables q,p,Si and S2: 

dq dH EO B 



dt dp ' 

dp _ _<9#eob 
dt dq 



= ^sT xSl ' "^sr xS2 - (l05) 

The knowledge of the time evolution of q(t), p(t), Si(t), Sa(t) is then injected 
into some gravitational-wave generation formalism, hij{t) = Hij(q,p, Si, S2), 



32 The PN-expanded EOB building blocks A(R), B(R), . . . already represent a resummation 
of the PN dynamics in the sense that they have compactified the many terms of the original 
PN-expanded Hamiltonian within a very concise format. But one should not refrain to further 
resum the EOB building blocks themselves, if this is physically motivated. 
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and used to compute the waveform h™ sp (t) during inspiral and plunge, up to 
some 'matching' time t m . [The analytical prediction Hij(q,p, S\, S2) can be 
computed with various accuracies: Ref. [3] used the lowest-order quadrupole 
approximation, while Ref. |120j uses a 3PN-accurate, resummed quadrupolar 
waveform.] Then, one continues this waveform to the merger and ring-down 
phase by smoothly matching, around the matching time t m , the EOB-dynamics 
derived /i^ sp (i) to a pure ring-down waveform, made of the superposition of 
several quasi-normal-mode frequencies: h lms (t) = 'J2c n e~ c ' n ( t ~ tm \ with a n — 

n 

a n + iuj n . Finally, this procedure defines a complete, quasi-analytical EOB- 
based waveform (covering the full process from inspiral to ring-down) as: 

h EOB (t) = 6(t m -t)H l] (q(t),p(t),S 1 (t),S 2 (t)) 

+ 9(t - t m ) C atchod e~ ff " (t -' m) , (106) 

where 9(t) denotes Heaviside's step function. 

The lowest approximation to the complete EOB waveform was constructed 
in [5] . Since then, more accurate versions were constructed in [1311 1120) . These 
EOB-type waveforms have been compared to full, 3d numerical relativity wave- 
forms in [1171 1129j . When taking advantage of the flexibility available in the 
EOB approach, an excellent agreement is reached between the quasi-analytical 
EOB-based waveforms and the numerical relativity ones. 

8 Conclusions 

In conclusion, we hope to have exemplified the way compact binaries set theo- 
retical challenges to General Relativity. 

On the one hand, over the past thirty years, binary pulsars have stimulated 
a lively dialogue between Experiment and Theory. This dialogue has led to 
novel tests of General Relativity, which have confirmed, with high accuracy, 
some of the strong-field and radiative aspects of Einstein's theory. The recent 
discovery of a double binary pulsar has greatly increased the number of available 
strong-field tests of General Relativity. 

On the other hand, the forthcoming detection of gravitational-wave signals 
in large interferometers is currently stimulating both analytical and numerical 
investigations in inspiralling and coalescing binary black holes. For the mo- 
ment, this fosters a dialogue between numerical results and analytical methods. 
Hopefully, one will soon be able to compare the combined analytical-numerical 
predictions to real gravitational wave data. 
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